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Abstract 

We study the asymptotics of recurrence coefficients for monic orthogonal polynomials 7r„(z) with the 
quartic exponential weight exp \—N (^2^ + ^iz*)], where i e C and e N, — > 00. Our goals are: 
A) to describe the regions of different asymptotic behaviour (different genera) globally in i e C; B) to 
identify all the critical points, and; C) to study in details the asymptotics in a full neighborhood near 
of critical points (double scaling limit), including at and near the poles of Painleve I solutions y{v) that 
are known to provide the leading correction term in this limit. Our results are: A) We found global in 
t € C asymptotic of recurrence coefficients and of "square-norms" for the orthogonal polynomials 7r„ for 
different configurations of the contours of integration. Special code was developed to analyze all possible 
cases. B) In addition to the known critical point <o = ^i^i found new critical points ti — and 
t2 — \. C) We derived the leading order behavior of the recurrence coefficients (together with the error 
estimates) at and around the poles of y{v) near the critical points ioi ii in what we called the triple scaling 
limit. Wc proved that the recurrence coefficients have unbounded C'(A^^^)-size (in t) "spikes" near the 
poles of y{v) and calculated the "universal" shape of these spikes for different cases (depending on the 
critical point io,i and on the configuration of the contours of integration). The nonlinear steepest descent 
method for Riemann-Hilbert Problem (RHP) is the main technique used in the paper. We note that the 
RHP near the critical points is very similar to the RHP describing the semiclassical limit of the focusing 
NLS near the point of gradient catastrophe that the authors solved in [5]. Our approach is based on the 
technique developed in [5]. 
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1 Introduction and main results 

In this paper we consider monic polynomials 7r„(z), orthogonal with respect to the quartic exponential 
weight e~'^^''^-*\ where f{z,t) = ^z^ + \tz^, t e C and N e N. As z -> cx), the weight function is 
exponentially decaying in four sectors Sj of the opening 7r/4, centered around the rays D,j — {z : argz = 
_aijt _^ 1^ j _ 1^2,3,4. We consider the most general case when the polynomials 7r„(z) are 

integrated on the "cross" formed by the rays fij, where the rays fli,2 are oriented outwards (away from 
the origin) and the rays fls^i - inwards. The corresponding bilinear form is 

{P, Q),ue2,ss,,. =Y.Si I P{z)q{z)e~''f^^''Uz, f{z, t) \z^ + Uz^ (1-1) 

j = l •'^i 

where Qj are fixed complex numbers chosen to satisfy Qi + Q2 ~ Pa + (?4- Moreover, since multiplying 
all the Qj^s by a common nonzero constant does not affect the families orthogonal polynomials, these 
parameters are only defined modulo the action of the group C* and hence the orthogonal polynomials 
are naturally parametrized by points in CP^. 

Alternatively, the bilinear form (1-1) can be represented 

as 

3 

(p, q), ^Y.'^k p{z)q{z)e-''f(^-'Uz, (1-2) 
k=l ■'^t' 

vi = -gi , V2 = -92 - Qi , = -Qi (1-3) 

where voj, j = 1, 2, 3, are simple contours emanating from 00 
along and returning to 00 along [3] (Fig. 1). 

Then in the case a) we have 1^2 7^ (and we therefore can, 
and will, normalize it to be z/2 = 1) and the following three 
cases are possible: 

1. The "generic case": viv^^ ^ 0, vi 7^ 1 ^ v^, so that 
there are three contours vjj in (1-2); 

2. The "consecutive wedges": either vi or (but not 
both) is zero so that there are two adjacent contours 
Wj in (1-2); 

3. The "real axis": = Q and vi = \. 

The remaining case b): Qi + Q2 = 0; corresponds to = Oj 
so that the following three cases are possible: 





















LU3 







Figure 1: The contours of integration and 
the asymptotic directions. The contour 
•074 is homologically equal to —wi — ^2 — 
•073, and it is unnecessary for the definition 
of the pairing (1-2). 
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1. The "single wedge": vi — 0, so that there is only one 
contour zu^ in (1-2), for which we can, and will, set 1/3 = 1; 

2. The "opposite wedges, generic": 1^11^3 7^ and vi ^ 1^3; 

3. The "opposite wedges, symmetric": vi — i^^ 0. 

The orthogonality condition for the monic polynomials 7r„(z) can now be written as 

{TTn,z'')p ^hnSnk, /c = 0, 1, 2, • • • , n, V = {vi , 1^2 , V3) , (1-4) 

where the coefhcient h„ can also be written as h„ = (7r„, 7r„)i7 and hence is the equivalent of the "square 
norm" of 7r„ (but it is in general a complex number). The existence of orthogonal polynomials 7r„(z) is not 
a priori clear. However, if three consecutive monic polynomials exists, they are related by a three-term 
recurrence relation 

7r„+i = (z - /3„)7r„(z) - Q;„7r„_i(z), (1-5) 

where a„ = an{t,N), /3„ = f3n{t,N) are called recurrence coefficients. 

If the bilinear pairing is invariant under the map z i— > — z then it follows immediately that the 
orthogonal polynomials are even or odd according to their degree and thus /3„ = , Vn G N (for example 
in case (al) with coefficients vi = v^, and 1^2 — 0). Then the remaining recurrence coefficients a„ satisfy 

n 

a„[l -f t{an+i + Oin + a„-i)] = — , (1-6) 

which is known in literature as the string equation or the Freud equation [15]. We are interested in 
the asymptotic limit of q;„,/3„ as iV ^ 00 and = a; > is fixed and finite, so we will use notations 
an = an{x,t), Pn = I3n{x,t) instcad of a„ = a„(t,iV), /3„ = l3n{t,N). 

In the case of V2 — —1 and x = 1 (that is, n — N) and a fixed t G {^j2^ 0)) the asymptotics of q;„, /3„ 
was obtained in [10] as 

a„(l,0 = ^^^P^ + O(n-i), (1-7) 

and /3„ decaying exponentially as n — >■ 00. (To be more precise, Theorem 1.1 from [10] states that there 
exists some no = ^o(0i such that a„(l, t), /3„(1, t) exists for all n > ng and have the above mentioned 
asymptotics.) In the non symmetrical case, however, the recurrence coefficients /?„ are, generically, 
different from zero. Then, instead of (1-6), we have the general Freud system (we indicate how to derive 
it in Section. 3.1): 

= /3„ + 1 [(2/3„ + /?„+i)a„+i + {Pi + 2a„(l - (5„o))/3n + an/5„-i(l - (5„o)] , (1-8) 
— ^ an+t [a„Q;„_i -|- a^,(l - 5„o) + an+i^n -I- Plan + a„^„-i(/3,i + Pn-i)] ■ (1-9) 
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where Sij denotes the Kronecker's deha. Assuming that /?„ — > /3 and a„ — >■ a as n — >■ od, we obtain two 
leading order algebraic equations 

(3{l + 6ta + tl3^) = , a{l + 3ta + 3tp^) = X, (1-10) 

which have two solutions 

Vl + I2xt - 1 , , 

13^0, a^^ , (1-11) 

and 

r.9 1 \/l - 15xt - 1 , , 

/32^-6a--, a=- . (1-12 

The dependence on x is rather fictitious: indeed, looking at the pairing (1-1) one sees that if < a; 7^ 1 
then we can rescale z = t — x t and obtain the case N = n {x = 1) without loss of any generality; we 
shall assume this done throughout the paper, but still distinguish n and N because they play a slightly 
different role. 

We point out that the asymptotics of the recurrence coefficients for orthogonal polynomials with in- 
tegration on the real axis (and analytic continuation thereof in the complex i~-plane) satisfies identically 
/?„ = 0, and hence only the first solution (1-11) is relevant. However, in view of applications to combina- 
torics of maps, it is not clear what role, if any, the second solution (1-12) play. It could be, perhaps, of 
some relevance that the critical point ti — oi the second solution is actually closer to the origin than 
the critical point io = — of the "standard" (first) solution. 

The recurrence coefficients problem was studied in [10], following earlier work [14], for negative values 
of f e {—j^,0). For definiteness we shall assign arg(t) — tt, so that (refer to Fig. 1) the contour wi consists 
of the two rays arg(z) = n72 of arg(z) — j,^ and ru^ of arg(z) — (other determinations of 

arg(t3) would only reshuffle the contours around). To describe the results of [10] and to set the stage for 
our results, we need to recall that the general solution of the Painleve I (PI) equation 

y"iv) = 6y^iv)-v (1-13) 

is parametrized by two parameters in terms of a certain Riemann-Hilbert problem described in Section 
2. It is known that any solution to PI has infinitely many poles with a Laurent expansion of the form 

The Painleve property asserts that the only singularities that can occur are of this form, that is, the 
position of these poles depends on the chosen solution, and it is largely unknown, except for some 
asymptotic localization of the remote poles, see, for example, [16]. In the following theorem and henceforth 
we use notations a„(i) = a„(l,i), /3„(t) = /3„(1, t). 
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Theorem 1.1 ([10]) Let y^^^v) :— y{v;l — i':}) , y^^^^ (v) := y{v;h'i) (see Def. 2.1 below), be the two 
solution of Pl(l-13) Let K <Z C be a compact set that does not contain any of the poles of y^^\y^^\ 
Let t E C approach the critical value to = —j^ in such a way that 

(1-15) 

where v E K . Then, for large enough n, the recurrence coefficients a„(i),/3„(t) have asymptotics 

anit)^2-2hiiy^^\v)+y'^'^\v))N-r. +0{N~'i), 

/3„(i) = 2To35(y(°)(w)-y(iHi;))iV"5 +0(iV-t) (1_16) 

as n = N ^ +00, which is valid uniformly in K . Moreover, the O terms can be expanded into a full 
asymptotic expansion in powers of . 

The statement of Theorem 1.1 is an example of what is known as the double scaling limit near a critical 
point and it is obtained using the steepest descent analysis and a special "Painleve I parametrix" that 
was first introduced in [14]: it does not address, however, the asymptotics of the recurrence coefficients 
when V is at or close to (in a "triple" scaling sense to be specified later) a pole of either y^^\v) or y^'^'(w). 
Also, no information is available for the case t < —1/12, as well as for general complex values of t. Thus, 
the main results of this paper are: 

1. finding the global (in i e C) leading order behavior of the recurrence coefficients a„(t), /?„(t) and of 
the "square- norms" h„ for the orthogonal polynomials 7r„ for the cases to different configurations 
of the contours zuj as listed above; 

2. deriving new critical points = in the case 1/2 = (this case was not considered in [10]) and 
t2 — \', note that ti is closer to the origin i = than t^; 

3. deriving the leading order behavior of a„(t),/3„(t) (together with the error estimates) at and 
around the poles of the PI solutions j/*^"'''^^^ (w) near the critical points to,ti; we will see that 
Cin{t) , j3n{t) and h„ have unbounded "spikes" near the poles of y^^^'*^^-' (w) and study the shape of 
these spikes in certain cases. 

In regard to the point <i = (which, to the best of our knowledge, was not considered in the literature), 
we also find an asymptotics that is related to the same Painleve 1 equation in the theorem below. 



^The parameters that were indicated with a, /3 in [10] in our notations are: o = 1 — fi, /3 = 1/3. 
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Theorem 1.2 Let y{v) := y^^'^{v) = y{v\vi) (see Def. 2.1) he the solution of PI (1-13); let K C C be a 
compact set not containing any poles of y^^\ Let t depend on N so that 



N^St := TVS U 

15 



where v € K (see 6-8). Then, uniformly for v € K , we have 



3t2t5' 



(1-17) 



ifise 10 3 6^6 10 3 

= -1 + % y{v) + 0(7V-5) , /3„ = -ii - 1, y{v) + 0{N~-^) 



Ni 



exp 



97V 195A^ 



TVs 



St + e" 



.6i 
N 



tHi 



(1-18) 
(1 + 0(^^-5)), (1-19) 



where Hj = ^(y')^ -\- yv ~ 2y^ is the Hamiltonian of PI (2-5) and H[{v) = y{v). 



Theorem 1.2 is, of course, of the same nature as Theorem 1.1. It is clear, though, that in order to 
study the full neighborhood of a critical point tj, j — 0, 1, which, by analogy with the zero dispersion 
limit of the focusing Nonlinear Schrodinger equation (NLS) will be called a point of gradient catastrophe, 
one must separate the asymptotic analysis in two distinct regimes: 

• Away from the poles: the variable v is chosen within a fixed compact set that does not include 
any pole of the relevant solutions to PI; 

• Near the poles: the variable v undergoes its own scaling limit and approaches a given pole at a 
certain rate. 

Theorems 1.1, 1.2 are examples of the regime "away from the poles". To investigate the regime "near 
the poles" we must use a novel modification that we could call triple scaling. 



Theorem 1.3 Consider the setups as in Thm. 1.1 with vi ^ l — v:^ or Thm. 1.2, with the same notation 
the former case and y^^^ in the latter. Let Vp denote any chosen pole of y^^^ (which is, if 
in the first setup, not a pole of y^^^^ ). Let t approach to ~ o?^ = in such a way that it satisfies 
respectively 

1 Vp s 1 Vpe~^ s 

^+12 ="]vt3f2f "sTItV °' ^~ 15 ^~3f2W2iV' 

Then 

hi 1 



'^-^^^ = i-^2+0{N'h^^), (1-21) 

/3„(t)=ao-f ^- ^, (1-22) 

^' 2s{l-hos) + 0{N-iy ^ ' 
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where ao,6o (the limiting values of a{t),b{t) from Table 2) are given by oq = 0, &o = 's/S and Cq 
— 3j , bo = 2i in the cases t ^ to and t ^ ti respectively. The numbers h„ satisfy: 



n2^ exp 



3A^ iVswp 



V2! 



h„ = 7r(-l) exp 



9N 13 N^Vpe- 



13 . 



352s 



2i - - +0(N-r.s-^) 
s 



15' 



(1-23) 
(1-24) 



These formulte hold uniformly for bounded values of s as long as the indicated error terms remain in- 
finitesimal. In particular s can approach s — or s — ^ at any chosen rate 0{NP) with p € [0, \) (the 
case p = allowing any given fixed value s 0, ). 



As the reader notices, the asymptotics has a dramatically changed form and does not involve now 
any transcendental function. Note that the scale of the phenomenon in this case is 0{N^^) aromid the 
location of the image of the pole Vp (see (6-7) or (6-8)respectively) in the t-plane, whereas the scale at 

4 

which the transcendental nature of the asymptotic is shown is N~ = . To study this new phenomenon, it 
is convenient to set a triple scaling of the form 

t = .o + ^ + |, (1-25) 

where the value of parameter Ci corresponds to a particular pole of the Painleve I transcendent. Note 
that, according to (1-21), (1-22), the values of a„,/3„ are unbounded as s — )■ and s ^ bo (the latter 
is valid only for /3„). A quite different phenomenon occurs instead if we are in the setup of Theorem 
1.1 with the same triple scaling limit but with additional symmetry i/i ~ 1 — (the case excluded from 
Theorem 1.3). In this case the two functions j/*^"^ and y'^^-' are the same solution to PI (1-13) and a sort 
of cancellation in (1-21), (1-22) occurs. 



Theorem 1.4 Consider the setup of Thm. 1.1 with t approaching to and vi = 1 — v^. Let Vp be a pole 

of y{v) '.= y^^\v) ~ y'^'^\v) and let t vary so that 

, . 1 Vp s 

t = tp{s) - - J^ - 2f3f^f + Wm' 

where s = 0{N^p) with an arbitrary p G [Oj ^)- Then the following holds: 

9 9- s^ + 0{N~l) 



4 1 



0{N- 



/3„ = 0, 



7r\/8 2 



N 



exp 



3A^ w„ 



3i2t 



3-s 
1 + s 



(1-26) 



(1-27) 
(1-28) 
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The variable s may approach the points s — ±1 at some rate (a quadruple scaling^ as long as the 
corresponding error indicated in the formultB above terms are infinitesimal. 



Remark 1.1 Note that the values of an is unbounded in the vicinity of s — ±1 and h„ is unbounded 
in the vicinity of s = —1 (there is no information for h„ in the vicinity of s — +1). Let us denote 
the Hankel determinants of the moments by An(t,N) (see Remark 3.1) and use tp{s) as in (1-26): since 
Oin = ^n-iA„4-i deduce that A„(t(s),n) vanishes at s ~ ±1 (within our error estimates), while 



A2 

A„±i(i(s),n) vanish at s ^ {lj3} and s G {—3,-1} respectively. 

2 The Riemann— Hilbert problem for Painleve I 

Let the invertible matrix-function P = P(^,w) 
be analytic in each sector of the complex ^-plane 
shown on Fig. 2 and satisfy the multiplicative 
jump conditions along the oriented boundary of 
each sector with jump matrices shown on Fig. 2. 
The entries of the jump matrices satisfy 



1 + LOoOJl = -W_2, 
1 + UJqUJ^i = -UJ2, 
1 + UJ^2^-1 = ^1, 



(2-1) 



1 


0" 




■ 1 


W 1 


_ W-2e 


1 







1 



so that the jump matrices in Fig. 2 depend, in fact, 
only on 2 complex parameters (that uniquely de- 
fine a solution to PI). The matrix function P(^, w) 
is uniquely defined by the following RHP. 



Problem 2.1 (Painleve 1 RHP [16]) The ma- 
trix P{^;v;lu) is locally bounded, admits boundary here •& :— {}{£^;v) :— — w^i 
values on the rays shown in Fig. 2 and satisfies 




Figure 2: The jump matrices for the Painleve 1 RHP: 



P+ = P M, 



(2-2) 
(2-3) 



where the jump matrices M — M{^]v,lj) are the matrices indicated on the corresponding ray in Fig. 2, 
with LJ :— (cj_2, 1^-1, Wq, wi, a;2) satisfying (2-1). 
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For any fixed values of the parameters w^, Problem 2.1 admits a unique solution for generic values 
of v; there are isolated points in the w-plane where the solvability of the problem fails as stated. The 
piecewise analytic function 

^{^,v;uj)^P{tv;uj)e''^'' (2-4) 

solves a slightly different RHP with constant jumps on the same rays and thus solves an ODE. Direct 
computations using the ODE and formal algebraic manipulations of series along the lines of [13, 11, 12] 
show that ^ admits the following formal expansion 



1 -i 
1 i 



V2 

Hia, Hfl+ya2 jv' - 'lH]-2y') iy'^2iHiy ^ \ 

V4 2^ 24^2 4^2 J 

Hi:=]^{y'f + yv~2y\ 7? := «) = _ , as C ^ oo. (2-5) 

where y = y{v) solves the Painleve I equation (1-13). The matrix \E'(^, -u; cD) uniquely defines a solution 
y{v;u}) of PI (1-13), and viceversa. The family of solution we shall use consists of the choice wq = in 
(2-1). Then the constant jump matrix for depends only on one free parameter: we shall choose it to 
be wi, with Lu±2 = — 1 and = 1 — uj-i. 

Definition 2.1 The functions y{v\ x) are the solutions to Painleve 1, which are defined via (^, u, (—1, >f, 0, 1 — x, — 1)) 
in Problem 2.1. We shall abbreviate this notation by {S,,v\ x) . 

3 The RHP for recurrence coefficients 

It is well known ([7]) that the existence of the above-mentioned orthogonal polynomials 7r„(2;) is equivalent 
to the existence of the solution to the following RHP (3-1). More precisely, relation between the RHP 
(3-1) and the orthogonal polynomials 7r„(z) is given by the following proposition ([10]), which has the 
standard proof (see [7]). 

Proposition 3.1 Let 57 :— ^^"^=1 '^'^'^ define v : ^ C by i'{z) = Vj when z € Wj. Then the solution 
of the following RHP problem 

Y{z) is analytic in C \ 51 

r+(z) = r_(z)(^ J ),zeii, ^3_^^ 

r(z) = (l + 0(z-i))(^ 'I ^z^oc. 

exists (and it is unique) if and only if there exist a monic polynomial p(z) of degree n and a polynomial 
q{z) of degree < n — I such that 

(p(z),z'=^)^^,^,,^3 =0, forall A: = 0,l,2,---,n-l, (3-2) 
((7(z),z'')j,j^^2^^3 = 0, for all A; = 0, 1, 2, • • • , n - 2, and {q{z), z"'~'^)^^^^^,^^ ^ ~2TTi. (3-3) 
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In that case the solution to the RHP (3-1) is given by 



p{z) Cn[p{z)v{z)e~''f(^''^] 
q{z) Cn[q{z)v{z)e-^1^^-^')] 



e C \ where CnW\ = 



1 

2TTi 



Y{z) = 

is the Cauchy transform of (j){z). 

Remark 3.1 It follows immediately that the polynomials p^q in Proposition 3.1 coincide with 



(3-4) 



T^n{z) = — det 



Mo Ml 

Ml 



Mn-1 

1 Z 



■ ■ M^^— 1 Mti 
Mn+1 



M2n-2 M2n-1 



liz) 



-2iTT 

"a7 



det 



Mo Ml 
Ml 



Mn-2 

1 Z 



■ Mn-2 Mn-1 
M« 



M2n-4 M2n-3 
^n-1 



respectively, where pLj (z^ , 1),^^ ,^2.1/3 are i/ie moments and A„ := det ^^^^ j^ . /i is clear from 

these expressions hut it is also a well known fact [6] that the condition of existence of the n-th orthogonal 
polynomial pn is that A„ ^ 0; on the other hand it is known from [I4] that existence of pn is equivalent to 
the solvability of the RHP and hence the existence of the solution for the RHP problem (3-1) is equivalent 
to A„ ^ 0. This determinant is sometimes referred to as the "tan function" of the problem [4j. Note also 
that the "square-norms" of the polynomials are ratios of Hankel determinants 

A„+i 



(7r„(z), 7r„(z)) - = (7r„(z), z") - 



If 7r„+i(z), 7r„(z) , 7r„_i(z) are monic orthogonal polynomials then they satisfy 

7r„+i(z) = (z - /3„)7r„(z) - a„7r„_i(z) 



(3-6) 



(3-7) 



for certain recurrence coefhcients a„,/3„ [18, 6]. The following well known statements (see, for example, 
[14], [7], [10]) show the connection between the RHP (3-1), the orthogonal polynomials 7r„(z) and their 
recurrence coefficients. 



Proposition 3.2 Let r(")(z) denote the solution of the RHP (3-1). If we write 





r(")(z) = 1 + + + o(z-3) 



then 



z^ 



/?„ = 



z — > 00, 



(") 



(3-8) 



(3-9) 
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Proposition 3.3 Suppose the RHP (3-1) has solution y(")(z). Let (3-8) be its expansion at 00 and let 
oimPn be given by (3-9). If an 7^ 0, then the monic orthogonal polynomials 7r„+i(z), 7r„(z) , 7r„_i(z) 
exist and satisfy the three term recurrence relation (3-7). 
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3.1 String equation for q;„(x, t), Pn{x,t) 

The string equations, or Freud's equations, for the recurrence coefficients an, (3n are nonlinear difference 
equations. Assuming that the corresponding orthogonal polynomials exist, they can be obtained as 
follows. On one hand we have 

Z7r„(z) = 7r„+i(z) + /3„7r„(z) + a„7r„„i(z). (3-10) 
One can iterate (3-10) to find z'^7r„(z) for any k eN. On the other hand we have 

EE ^ I., /" (^„7r„,e^^/) dz = ^ i/, / (7r>™ + n^TT'^ - iV7r„^„/' (z)) e'^^dz = (3-11) 

= + {TTn,n'^)^ - 7V(7r„,/'(z)7r„)^ (3-12) 

Since f'{z) is a polynomial, the last term above can be written as a polynomial in the recurrence coeffi- 
cients using repeatedly (3-10). For n = m the first two terms are the same and vanish because tt^ is a 
polynomial of degree n—1 and 7r„ is orthogonal to any polynomial of lower degree. Then (3-12)„„ yields 
a recurrence relation. For m = n — 1 we have 

«, 7r„„i)^ + (7r„, - N (7r„, /'(z)7r„_i)^ ^ n (7r„_i, 7r„_i)^ - N (tt,,, /'(z)7r„_i)^ . (3-13) 

Equation (3-13) yields (1-9) while (3-12) for n = m yields (1-8). If the orthogonality pairing is symmetric 
under z i— )■ — z, that is, if 

{p{z),q{z))^^{p{~z),q{^z))^ (3-14) 
then it follows easily that /3„ = and then (1-8, 1-9) reduce simply to (1-6). 

4 Steepest descent analysis of the RHP (3-1) 

The steepest descent analysis in general terms for these kind of orthogonal polynomials with a polynomial 
external field was investigated in [3] and so we refer the reader there for details. The schematic of the 
approach is outlined here; as customary, the problem undergoes a sequence of modifications into equivalent 
RHPs until it can be effectively solved in approximate form while keeping the error terms under control. 

• One starts with the problem for Y (3-1) and seeks an auxiliary scalar function g{z), called the g- 
function, which is analytic except for a collection S of appropriate contours to be described subsequently 
and behaves like Inz -1- 0{z~^) near z — oo: the contour fl of the RHP (3-1) can be deformed because 
the RHP (3-1) has an analytic jump matrix. The final configuration of Q must contain all the contours 
where g is not analytic. 

• Then we introduce a new matrix 

T(z) e-^^^F(z)e-^(s(^^*)-i)-^ (4-1) 
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As a result, T solves a new RHP 
T{z) is analytic in C \ 

< r+(z)=T_(z)l ^ ef(H-/.-) j.^ef}, where ;i(z,t) :=2g(z,t)-/(z,t)-l 

r(z) = (l + 0(z-i)), z^oo. 

(4-2) 

• At this point the Deift-Zhou method can proceed provided that the function g{z), the constant £ and 
the collection of contours into which we have deformed the problem fulfill a rather long collection of 
equalities and -most importantly- inequalities that we set out to briefly describe [20, 3]: we say here 
that if all these requirements are fulfilled the full asymptotic for the problem can be obtained in terms 
of Riemann Theta functions on a suitable (hyper) elliptic Riemann surface of a positive genus (with the 
case of zero genus not requiring any special function). 

4.1 Requirements on the ^f— function 

The (deformed) contour Q can be partitioned into two disjoint subsets of oriented arcs that we shall 
denote by DJl and term main arcs, and £ or complementary arcs; this partitioning is subordinated to 
a list of requirements for g and h. 

4.1.1 Equality requirements for g 

1. g{z) (to shorten notations, we drop the t variable in this subsection) is analytic in C \ (9Jl U €.) and 
has the asymptotic behaviour 

g{z) = In z + 0{z-'^) , z ^ oo; (4-3) 

2. g is analytic along all the unbounded complementary arcs except for exactly one unbounded com- 
plementary arc which we will denote by 70, where 

g+ (z) - .g_ (z) = 2i7r z e 70 (4-4) 

(note that the function e^^*^^) is analytic across all the unbounded complementary arcs since N = 
n G N by definition); 

3. on the bounded complementary arcs 7c, the function g has a jump 

5+(z) - .g-(z) = 2Trir]c , 77c G K , z e jc, (4-5) 
where rjc is a constant on each connected component of the complementary arcs; 
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4. across each main arc (which are aU bounded by assumption) we have the jump 

zem. (4-6) 

We stress that the constant £ is the same for all the main arcs. 

Assuming that the contours dJl, £ are known, the function g{z) can be considered as the solution of 
the scalar RHP, defined by conditions 1-4. Similarly, the function h = 2g ~ f can be considered as the 
solution of the scalar RHP with the jumps 

^{h+{z) ~ h_{z)) ^ 2TTir]c, z e VJl, ^ {h+{z) + h_{z)) = 0, z G £, h+{z) - h_{z) = Ain, z E Jo, 

(4-7) 

and the asymptotic behavior 

/i(z) = -/(z) + 2 In 2 + 0(2-1) , z^oo. (4-8) 

It follows immediately from (4-7) that 5R/i(z) is continuous across the complementary arcs C 

4.1.2 Inequality (sign) requirements (or sign distribution requirements) for h and the 
modulation equation 

1. along each complementary arc 7c wc have 3?/i(z) < 0;with the equality holding at most at a finite 
number of points. In the generic situation these would be only the endpoints (we shall call this case 
regular, with the same connotation as in [8]); 

2. on both sides in close proximity of each main arc 7™ C StJt we have 3?ft.(z) > 

The sign distribution requirement for the main arcs implies that dih{z) is continuous everywhere in 
C and the main arcs belong necessarily to its zero level set. The main arcs 7™ can be considered as 
the branch-cuts of a hyperelliptic Riemann surface 5K(i), associated with g and h. The number of main 
arc (the genus of fH(t) plus one) needs to be chosen in such a way that the above sign conditions will 
be satisfied. The location of the endpoints A of each main arc (which are the branch-points of 9^(t)) is 
governed by the requirement 

m{z) ^0{z- X)^) , z^A, (4-9) 

known as the modulation equations. Since the jumps on the complementary arcs are constants, the above 
requirement can also be stated as 

h'{z) ^ OiVz - A) , z A, (4-10) 

where the discontinuity is placed on the main arc. The logic behind all the above requirements and the 
modulation equations will be briefly discussed in Subsections 4.2, 4.3. Note that the modulation equation 
(4-9) implies that there are three zero level curves of 3?/i emanating from each branch-point A. 
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4.1.3 The g-function and the modulation equations in the genus zero case 

Due to the modulation equations 4-9, solutions of tlie RHPs for g and for h commute with differentiation. 
Thus, the (scalar) RHP for g' is: 

1. g'{z) is analytic (in z) in C \ 9Jt and 

5'(z) = - +0(z-2) as z^oo; (4-11) 
z 

2. g'{z) satisfies the jump condition 

g'+ + 9- = f on m. (4-12) 

Let us consider the case of a single main arc 7^ with the endpoints Aq, Ai. Using the analyticity of /', 
the solution of the latter RHP is given by 

.'(.) = ^ (c/.]g(C)/ ^ ' ^(^^ - , (4-13) 

where the contour 7^ encircles the contour 7^ and has counterclockwise orientation (z is outside 7m). 
It is known ([10]) that the case t & {—^^,0) is the genus zero case with real branch-points (we will derive 
the same result shortly). Using (4-13), the asymptotics (4-11) yields two equations 

called moment conditions, which are equivalent to the endpoint condition (modulation equation) (4-9). 
We use the moment conditions (4-14) to define the location of the endpoints Aq.i, where we put Aq < Ai. 

In the case of a polynomial f{z,t), equations (4-14) can be solved using the residue theorem. Setting 
Aq = a — 6, Xi ~ a + b and using the residue theorem on equations (4-14) we obtain 

a + ta{a^ +'^b^) ^0 and + ^b^ + t{a* + 3a^b^ + "^b^) = 2. (4-15) 
2 2 8 

There are two possibilities: a = and a 7^ 0. In the first case we obtain solutions to the system (4-15) as 

a-0 52 = -^(lT%/rTl2t), (4-16) 

Aoa = Tb = T^-^^{l-VTTm) (4-17) 

The choice of the negative sign in (4-16) comes from the requirement that b is bounded as i — > 0. Observe 
that for t > to — — j;^ the values ±6 coincide with the branch-points, derived in [10]. At the critical point 

i = io = -^ (4-18) 
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the two pairs of roots (4-16) coincide, creating five zero level curves oi^h{z) emanating from the endpoints 



±60, where 60 = V8- The second pair of roots b = ±y^-§(l + VI + 120 are sliding along the real axis 
from ±00 to ±60 as real t varies from 0^ to to, and sliding along the imaginary axis from ±ioo to as 
real t varies from 0"'' to +00. 

In the second case a ^ 0, the system of modulation equations (4-15) becomes 



t{a^ + ffo^) = -1, 
t62(a2 _ |52) = 2, 



which yields 



= --^(1± Vl - 150 and ^ -^(3^2^1-150- 



(4-19) 
(4-20) 



4.1.4 Explicit computation of g and h. 

Once the values of branch-points (endpoints) Ao,i are determined, one can calculate explicitly g(z) and 
h{z) = h(z;t), where h{z) — 2g{z) — f(z) — I. The expression 

Ri^) f no 



h'{z) 



2m 



(4-21) 



iC-z)R{0+ 

for h'(z) is readily available from (4-13) by placing z inside the loop 7™- However, it seems easier to 
calculate h'{z) explicitly by solving the scalar RHP that h'{z) satisfies: 



1. h'{z) is analytic (in z) in 7^ and 

h'{z) ^ -z-tz^ + 2z^^ + 0{z^'^) as z ^ 00; 

2. h'{z) satisfies the jump condition 

h'_^ + h'_ =0 on 7m, 



(4-22) 



(4-23) 



which can be easily obtained from the RHP (4-12) for g'{z). There are two cases, symmetric and 
nonsymmetric depending on the value a ^ or a =/= 0. 

Symmetric case: a = 0. The RHP for h' has a unique solution (with h'^. £ L'^{jm)) that is given 
by h'(z) = —{k + tz'^)^/ z'^ — , where the endpoints ±6 of 7,„ are known and the constant k is to be 
determined. Assuming that b^ is given by (4-16), we obtain fc = 1 + ^tfe^, so that 



h'{z) = 



tz^ + 1 



{z^-b'Y-=- 



tz' 



Vl + 12< 2 



2^ 2 



(4-24) 



Since the branch-cut of the radical is [—6, b] we conclude that h'{z) is an odd function. Direct calculation 
yield 

h{z) = 2 In ^^ + ^f' _ ^ (2^2 ^ ^^2 + 4) (^2 _ i ^ (4_25) 
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It is clear that h{b) = 0. There is the oriented branch-cut of h{z) along the ray (— cx),— 6), where 
h+{z) — h^{z) = AiTi. Combined with (4-24), that implies 

h{z) = 0{z-h)i (4-26) 

(or a higher power of (z — h)). At the point of gradient catastrophe io = ^ j^, the order 0{z — b)i in 
(4-26) should be replaced by 0{z — 6) S . Because of the 47ri jump along (— oo. — fe), the function h{z) does 
not have 0{z + b)i behavior near z — —b; however, ^h{z) does have 0{z + b)^ behavior near z — —b. 

Non-symmetric case: a 7^ 0. Following the same lines (the algebraic computation being a bit more 
involved) one obtains 

h'{z) = -t (z^ + az- b^) ^{z-aY-b^ (4-27) 
where a, b are given by (4-20). A direct computation yields va w — z — a as (using tb^a^ — ^tb'^ ~ 2) 



h{z) ^ / h'{C)SC = -Vw^-b^ 



-b 



t 2 

-{w'^ - b'^){w + Aa) + j^w 



21n( " + ^f ) . (4-28) 



Remark 4.1 One can verify directly that h{z) satisfies the following RHP: 
1. h{z) is analytic (in z) in C \ {'y„i U (—00, Ao)} and 



h{z) ^ -^tz* -^z'^ + 2\nz- £ + 0{z-^) as z ^ 00, (4-29) 



whe 



9 9 1 

£ — \n — (symmetric case, with b given by (4-16).) (4-30) 

6^ 2a2 + 52 1 

£ = In , (nonsymmetric case, with a,b given by (4-20) (4-31) 

2. h{z) satisfies the jump condition 

h+ -\- = on 7m and ~ — Airi on (— oo,Ao), (4-32) 

Remark 4.2 As it was mentioned above, the solution to the scalar RHP for h commutes with differenti- 
ation in z; on the same basis, it commutes with differentiation in t as well. Thus, we obtain the following 
RHP for hf (symmetrical case): 

1. ht{z) is analytic (in z) in C\ 7™ and 

ht{z) = -^z"^ - it + 0{z-^) as z^oo, (4-33) 



where 



it = = 2 + 12t- 271 + 12^ . 3 
* 32 24<2 ' ^^■^^> 
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2. ht{z) satisfies the jump condition 



ht+ + ht- =0 on 7„. (4-35) 

These RHP has the unique solution 

ht (z)^-^ (2z2 + 62)^^2 _ ^2 (4_36) 
8 

that can be verified directly. In the non- symmetrical case, ht can be calculated in a similar way. 
Remark 4.3 The g-function gt{z) was defined in [10], eq. (3.2), as 

gt{z)^ I Hz~^)d^lt{0, (4-37) 

where Ht is the equilibrium measure in the external field f{z,t) and x — n/N = 1. In the case t >0, the 
equilibrium measure fit is the unique Borel probability measure on M that minimizes the functional 

If(p)^ J J lnJ^^^d^lix)d^l{y)+ J f{x,t))d^iix) (4-38) 

among all Borel probability measures on M. (Here the subindex t does not mean differentiation. ) The 
measure fit can be calculated explicitly. It turns out to be supported on the interval I = [—6,6], where 6 
is defined by (4-16), and it has a density given by 



dfi _ t f ^ Vl + 12i + 2 
dx 2tt 



(^x^ + "^^^af V^^, X G hb, 6]. (4-39) 



(In fact, the case t G (io,0), the equilibrium measure fit minimizes (4-38) among all Borel probability 
measures with the support on [—6, 6] .J The function gt{z) satisfies the requirements of Section 4-1-1 
because of (4-37). Since the RHP for g has a unique solution, we have gt{z) — g{z). 

4.2 Discussion about existence of h 

To the reader it could be a little bit of a mystery as to why there exists any function h{z,t) satisfying 
the above long list of conditions. However, this result was proven in a general setting, that is, for any 
polynomial f{z,t) and for any t G C in [2]. The idea of the proof is quite simple. Suppose that we 
have our contours fl and we want to find the h{z,t) function for a specific value t G C. Assume, on the 
other hand, that for a certain value of the parameter t (for example, for G {—j2i 0),) one can somehow 
find h{z,tsf), satisfying all of the above requirements (for example, h can be calculated directly using the 
residue theory, as above, or by use of the potential theory). Then one chooses a path in the parameter 
space (t-plane) that connects to t and shows that the requirements can be maintained throughout the 
path; we shall call this the continuation principle in the parameter space. This idea (implemented 
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in slightly different form) was at the basis of the discussion of [20] and [2]. In a general situation with 
f{z) being an arbitrary polynomial, the existence of a suitable h{z,t^) was established in [2], but in the 
present paper we will prove all the inequalities for h{z, t^), e { — j^,0), directly. In fact, the continuation 
principle is not limited to the polynomial or even rational potentials /. For example, in the context of 
the semiclassical limit of the focusing NLS, the continuation principle for a large class of analytic /, was 
stated and proven in [19] 

To indicate the obstacles that make the continuation principle nontrivial we point out that, as we 
follow our path in t-plane, it may happen that the regions where — < (the sea) moves in such a way 
to either pinch off one of the complementary arcs or to "expose" one of the main arcs (or causeways) ; 
in that case we can use local analysis to guarantee that a new main arc (causeway) or complementary 
arc respectively can be "sewn in" in order to adjust the situation; such an adjustments increases the 
genus of the solution. In fact it is quite a daunting task to try and describe in words this process; we 
invite the reader to have a close look at the pictures of the "phase diagrams" Figs. 5, 6, 7, 8, 9, 10. The 
reader should try and imagine how the main arcs and complementary arcs (which are not marked in the 
pictures) deform as we cross the phase-transition curves, also known as breaking curves, indicated there. 
In fact an interactive exploration tool was designed in Matlab and it is available upon request. 



4.3 Schematic conclusion of the steepest descent analysis 

The final steps in the steepest descent analysis involve adding additional contours, the lenses, which 
enclose each main arc and lie entirely within the — < region (the sea). One then re-defines T{z) 
within the regions between the main arc and its corresponding lens by using the factorization 



a d 




1 







of the jump matrices of T so that 

r+(z) = T_(z) 







1 

ad-^ 1 



i/-ie-^''+ 1 



(4-40) 



(4-41) 



where Vm is the (constant!) value of viz) on the main arc under consideration. Therefore, defining T{z) 
as T outside of the lenses and by 



T{z) ■.= T{z) 



(4-42) 



in the regions within the lenses and adjacent to the ± sides of 7^ one achieves a new problem with jumps 
that are constant on the main arcs and exponentially close to the identity or constant jumps on the lenses 
and complementary arcs. 

We spend a few more words for the "genus zero" case, namely, when there is a single main arc 7^ 
connecting two endpoints Aq, Ai, since this is the situation mostly relevant to the analysis here; the case 
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with several arcs, for the case of real potentials on the real line was fully treated in [8] and in the complex 
plane in [3]; while not being conceptually more difficult, it requires the introduction and use of special 
functions called Theta functions. 

4.3.1 The "genus zero" case 

This is the case when there is a single main arc "f„i that connects two endpoints Aq, Ai ; since the coefficients 
i^j are defined up to multiplicative constant, we can and will assume without loss of generality that they 
have been normalized so that the iy„i on the main arc satisfies Vrn = 1- Then the RHP for T is 

T-|_(z) = T^{z) ( -^Nh ? ) on the upper and lower lips respectively, 



r+(z)=T_(z)( l]=T_{z)t<j2 on 7^. (4-43) 



1 

-1 

r+(z) =T_(z) ( 1 ) on f]\7„ 

Due to the sign requirements, the off-diagonal entries of the jumps on the lenses and complementary axis 
tend to zero exponentially fast in any L^'-space of the respective arcs, p > 1, but not in L°° because at 
the endpoints Ao,Ai we necessarily have SR/i = 0. Near these points one has to construct explicit local 
solutions of the RHP called parametrices [8]. The type of local RHP depends on the behavior of h{z) 
near the endpoints. 

In a generic situation one has e^''^^^ = qNC^„'> {t)i(z~\j)2 {i+o{z~\j)) ^ j = 0,1 with cl^^\t) some nonzero 
constant. The critical case (or "gradient catastrophe" case) correspond to those special case whereby 
C^j)(t) =0 at one or the other or both endpoints, and thus 

where C['\Q is now nonzero (nondegenerate gradient catastrophe). In the former case the local 
parametrix can be constructed in terms of Airy functions and its construction is very well known since 
[8] (see also [10], [3]). The latter case requires the solution of a special RHP which can be reduced to an 
instance of the RHP for the Painleve I Problem 2.1. This was done in [10] and will not be repeated here. 
We point out that one of the main distinctive features is that 

• in the generic case there are three level curves SRft. = that emanate from the corresponding 
endpoint Xj (one of them being the main arc), see Fig. 3; 

• in the critical case there are five level curves = 0, one of them being the main arc (see Fig. 3) 

The final steps in the approximation mandates that we fix two disks Do,Di (small enough not to 
enclose any other endpoint) around the endpoints Aq, Ai and define a suitable approximate solution 

^ext{z) for z outside Dp.i 
$(z):-<( <i>eMMz) inside Do ' (4-45) 

<S>ea:tiz)riiz) insidc Di 
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Figure 3: The two typical configurations of level curves and sign distributions near the endpoint in the 
generic case (left) and critical case (right). Indicated are the complementary arcs 7c (there might be only 
one complementary arc in the critical case, depending on the function v^z)) and the lenses. The blue 
(darker) color corresponds to the region where — < (the sea). The parameter >c equals i^i or 1 — 
depending on the endpoint under consideration (see Fig. 11 and parameters therein). 

such that the error matrix £{z) :— T(z)$^^(z) solves a small-norm Riemann-Hilbert problem (as N 
oo) and thus can be -in principle- be completely solved in Neumann series. Here by 7^0,1(2:) we denote 
the parametrices near the endpoints Ao,i respectively. 

In all situations the matrix ^ext ("model solution" or "exterior parametrix") solves a RHP of the 
form (model problem) 

^extiz)+ = 'Pextiz)-i(T2 , 2 € 7m = [^0 , Al] , , ^ext{z) = 1 + Oiz~^) , Z ^ OO (4-46) 

with some particular growth behavior near the endpoints which depend on the scaling limit under con- 
sideration. In the usual case it satisfies 

^ext{z)^0{z-\j)-'i , z^A,, (4-47) 

but in special cases the behavior needs to be modified. 

At any rate, once we have achieved a suitable approximation for T{z), the recurrence coefficients for 
the orthogonal polynomials can and will be recovered via the formulae 

h„ = -2z7re^^ (ri)i2 , = iTi\, {T,\, , /?„ = - {T^\^ (4-48) 

where T{z) near 00 equals T{z) (since we are in the exterior region) and has expansion 

f (z) = T{z) = 1 + — + ^ + ..., z^oo . (4-49) 
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The latter coefficient matrices can be obtained from the corresponding expansion of ^extiz) near infinity, 
to within the error determined by £; in the generic (regular) case (t ^ tQ,ti,t2 and not on the breaking 
curves) the parametrices Vq, Vi are the well-known Airy parametrices and the standard error analysis 
(which we do not report here) shows that £ introduces an error of order 0{N~^). 

In this case the exterior parametrix (model solution) ^ext in the genus region is the "standard" 
solution (that we shall denote by 'So) to the following "model RHP": 



5*0(2) is analytic in C \ [Aq, Ai], 

*o+(^) = *o-(2)icr2 on [Ao,Ai], 

'I'o(z) = 1 + 0(z"i) as z 00, 



(4-50) 



The solution to the RHP (4-50) is given by 



*o(.) = ^^f^V(-3 + ..)^f^) , (4-51) 



2 V^--^o, 

which has expansion (recall our notation Ai = a + 6, Aq = a — 6) 

b abao 



*o(z) = l--^2 + ^l-^+0(^-') , z^^. (4-52) 
Thus, near z — 00, one finds 

T{z) - (^1 + io(iV-i)^ 'f.iz) ^ T, = (1 + 0{N-^mo.j, (4-53) 
where denote the Taylor coefficients of (z) at infinity. 
4.3.2 Recurrence coefficients in the genus cases 

As explained in Section 4.1.3 there are two types of genus zero solutions and hence the final formulas 

are different. Using (4-48), the approximation (4-53), the explicit form of V^o (4-51) and the explicitly 
calculated expressions for Ai, A2, one finds the results summarized in Table 1. 



.3.3 The regions of higher genera 



From the global analysis reported in Figures 5, 6, 7, 8, 9, 10, the reader can see that there are regions 
where the hyperelliptic surface of h'(z) has genus 1 or 2. In this case, while the general scheme of 
the steepest descent analysis remains intact, the solution of the relevant model problem for ^0 requires 
Riemann-Theta functions. Formulae can be found in [8, 3]. The recurrence coefficients also are expressible 
in terms of Theta functions. In fact the formulas in [8] could be directly applied here, simply by modifying 
the choice of the a, ^-cycles (in the standard lore of Riemann surfaces) as described extensively in [3]. We 
shall not write explicit formulae here since it would require setting up a good deal of additional notation. 
Suffice it to say that the nature of the resulting expressions is one of rapidly oscillating functions of N 
and t, with amplitude that depends only on t. 
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Symmetric genus {St := t + j^) 



Non symmetric genus {5t t ~ j^) 



h„ = 2.(^^^jexp 



N 



4- {Vl2St+iy 
2M 



h„ = 2vr^^ exp 



N 



9 + AiVl5St - 30St 
60t 



&2 _ y/l25i - 1 



62 _ iVlM - 1 



/3n = a = 



a= -i\ —(3 + 2iVl5Stj 



Table 1: The leading order approximations of the "square- norms" and recurrence coefficients in the two 
genus-zero cases: all expressions are understood to within an error term of 0{N^^). T he expressions for 
h„ = nbe^^ and ^'s are in (4-30, 4-31). In fact a more careful analysis shows that /?„ in the symmetric 
case is exponentially small [10]. The reason is that the RHP can be seen to be close exponentially to a 
RHP with a symmetry z 1— >■ — z, for which the expression for /3„ automatically yields zero. Note that there 
are two choices of signs for a in (4-20) (the choice of signs for 6 amounts only in exchanging the labels 
of the branch-points) that lead to different (but quite similar) formulae and results; we will formulate all 
the results for this particular choice whereby a ~ —Si. 



4.4 Contour deformation. 

A general discussion of contour deformations once the appropriate (/-function has been found can be read 
in [3] and [2] ; we give here a brief sketch. We advise the reader to accompany this part with the pictures 
that are provided plentiful. 

In general, the contours of integration for the pairing (p, q) p can be deformed by use of the Cauchy 
theorem: any deformation that we shall allow must be such that the deformed contour approaches 00 
along the same direction arg(z) = — | arg(i) + of the original contour, so as to preserve integrability 
(we may even mandate that each contour is a straight line outside of a sufficiently large circle). Indeed, 
from ?R.h{z) = —di {f{z) — 2g{z) + £) and from the fact that g{z) is bounded by a logarithm, we see that 
for \z\ large enough the sign of 5ft/i is the same as of — for which the above directions are the directions 
of the steepest descent. 

The final deformation of the contours must fulfill the following requirements, that we describe referring 
to the regions — 5ft/i > as (dry) land, — < the sea (or other watery expression) and the main arcs 
(where 3?/i = 0) as bridges or causeways: 

• Along each contour —dih is always nonnegative, — Sfi/i > 0, i.e. each contour does not get wet; 

• if two or more (oriented) contours have been deformed so that they go through the same bridge 
(main arc), then the traffic (i.e. the weight of that part of contour) is the (signed) sum of all the 
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traffics. For example if tui, zu2 are deformed so that they pass trough the same main arc, then the 
weight of that arc shall be vi — 1^2] 

• each bridge (main arc) must carry a nonzero traffic, or else one needs to find a different g-function; 

• the precise form of the deformed contours as they enter/leave a bridge (i.e. the complementary arcs) 
is largely irrelevant, but for definiteness we shall stipulate that they proceed for a short distance 
along the steepest ascent line or —dih. 

In order to offer some rigorous study we consider in more detail the symmetric case of genus 0. 

Lemma 4.1 In the case V2 = —1 and t G (— ^^^,0), the function ^Rh{z), where h{z) is given by (4-25), 
satisfies the sign conditions along the contour fJ. 

Proof. First, it follows from (4-25) that Rh{z) = on [-6,6]. To show that 5R/i(z) > immediately 
above the main arc [—6,6], it is sufficient, by the Cauchy-Riemann equations, to show that 3/i'(z) < 
on the upper shore of [—6,6]. The latter follows directly from (4-24). We can now use the oddness of 
h' [z) to show that 3ih{z) > also below the main arc. So, the correct signs around 7^ are proven. The 
correct distribution of signs of SR/i along the complementary arcs follows from the topology of zero level 
curves of 5ft/i. Because 'Rh{z) is even ( \z + z"^ — h'^\ is even), it is sufficient to consider level curves 
only in the right halfpane. Direct check shows that both terms in (4-25) have positive real part on iM"*". 
There are two legs of zero level curves of 3?/i emanating from z — b and four legs coming from infinity 
with asymptotes ±|, see (4-29). Denote these legs as x±j: i = 1,2 respectively. Since Rh{z) < 

as real 2; — > 6 -I- and Rh{z) > as real z +00, we conclude that Rh{z^,) = at some e (6, +00). 
So, the only possible topology of the level curves x±j is that xi is connected with x-i through z^ and 
X2, X-2 are connected to b (since Rh{z) is a harmonic function, its level curves cannot form bounded 
loops). Thus, one can choose as complementary arcs any smooth curves "on the land" between xi and 
X2 and between x-i and X-2, i-e., in the region where Rh{z) < that connect 6 and 00. Fig. 5 shows 
(in red) main arcs 7™, but not complementary arcs jc- However, level curves Xj can be visualized in the 
"snapshots" of 2-plane that correspond to t e (— j^,0). Q.E.D. 

Remark 4.4 Similarly to Lemma it is easy to establish the correct sign distribution outside jm 'in 
the case when h is given by (4-28), which is valid, for example, when 1^1 = 1^2 = 0, 1^3 =/= (Single wedge) 
and t G (0, ji^), see Fig. 7. For t G (0, both b^ and a^ are negative, so that a,b E iC. From (4-28) it 
follows that Rh{a) — and setting the orientation of upwards, we see that Rh'j.{a) ^ respectively. 
Thus, we have the correct sign distribution outside 7™. 
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5 Breaking curves and global phase portraits 



A breaking curve A in the complex i-plane separates the regions of different genera in the asymptotic 
behavior of the recurrence coefficients, or regions of the same genus but with different number of main 
arcs (see, for example, the breaking curve that joins t = —j^tot—jm Fig. 5, which separates two 
regions of genus 2). It satisfies the system of equations 



which is the system of 3 real equations for two complex variables z and t. Here the subindex k in hk 
indicate the genus of the Riemann surface £H(i) where hk is defined. In our cases the genus can be 0, 1, 2. 
To simplify notations, we will drop the subidnex k whenever the genus of h is obvious. 

We will consider the breaking curves in the t plane where the sign requirements fail because a saddle 
point z of (a point satisfying hz{z,t)—0) collided with the contour = 9Jt U £. That means that 
either a complementary arc is pinched by the rising "sea" or a main arc (causeway) is touched by the 
dry land because of the receding "sea". In any case, equations (5-1) will be satisfied at z = z. There are 
three cases of breaking that we consider: 

• genus symmetric, i.e., /io(z) is given by (4-25); 

• genus non-symmetric, i.e., hQ{z) is given by (4-28); 

• genus 1 (symmetric). 

The resulting equations when plugging the expression for h into the system (5-1) are relatively simple 
and could be analyzed analytically; we find it much more effective and informative to study and plot 
them numerically. 

5.1 Genus symmetric 

Using (4-24), we obtain the following equation for the saddle point: 



Substituting this into ^h{z) = and using (4-25) after some algebra yields the following implicit equation 
for the breaking curve 



h'l^iz) = and ^hkiz) = 0, 



(5-1) 




or 




(5-2) 



In 



1 + 2V1 + 12t + y^3[(\/l + 12t + 1)2-1] 



+ 



3[(Vl + 12t + 1)2-1] 



= 



(5-3) 



1 - Vl + 12i 



m 
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or 



In 



= 0, (5-4) 



where u — ^/l + 12t. Note that according to (5-4) A is a bounded curve that starts at the point of 
gradient catastrophe u — because for large u the expression in (5-4) tends to ln(2 -I- VS). To obtain the 
asymptotics of the breaking curve A near t — to, we use expansion (6-2) of h{z,t) near the branch-point 
z = b, where the coefficients Cq, Ci are given in Table 2, left column, to write 

h{z)^-^{z'-b')i + -^{z'-b')i+--- . (5-5) 

{z + b)2 {z + b)2 

This expansion is a direct consequence of the modulation equation. Using (5-2), we calculate — 6^ = 
^"^J^^"^^' _ y/i'^it^ to) ^ Since for t close to to the saddle point z{t) (that satisfies hz{z{t),t) = 0) is close to 
6, we have 

M.(^),t)^ ^°[^^(^'e]^^^[^^-(^-^)]^ o((t-to)l)^ [^^(^-^°)]! (iot5^+4)+o((^^Mi), 

^^'^ {-2bt)i {-2bt)i 15(25)3(-i)i^ ' 

(5-6) 

where we utilized the formulae for Cq, Ci. Now the requirement "^h = yields | arg(i — t^) — + Trk, 
/c S Z, so that the breaking curve A near the point of gradient catastrophe is tangential to 

27r 47r 

&Yg{t^to) = ±— + —k. (5-7) 

Note that there are various branches to keep track of: the principal branch of the radical ^Jl + 12t leads 
to the curve joining t~—-^iot~Q{u^Oiou = \ correspondingly), light curve from — to on Fig. 
4; the secondary branch leads to the curve that joins t — — j^^tot— \ {u — Qiou — 2 correspondingly), 
on Fig. 4. 

5.2 Genus non-symmetric 

In this case we are looking for zeroes of h'{z) satisfying + az — 6^ = 0, see (4-27). They are given by 



1 I /b^^ 



Substituting (5-8) in 3fJ/i(z) = 0, where h is defined by (4-28), and repeating the previous arguments, we 
obtain an implicit equation for the additional breaking curves. Leaving the lengthy but straightforward 
details aside, we obtain the curves on Fig. 4 that join t = to t ^ and t = IE respectively. 
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5.3 Genus 1 symmetric 

For the case of genus 1 there are 4 branch-points; the only situation where we can have the saddle point 
h'{z) = on the zero- level set is when the saddle point is between to distinct connected components 
of the zero level-set of 3?/i(z) = 0. It is seen from the modulation equations that /i'(z)^ is always a 
polynomial of degree 6; we look here for solutions where (/i'(z))^ is an even polynomial. Since we are 
seeking a solution of genus 1, there must be a single double root. By the symmetry this root must be at 
the origin; this allows us to write 



h'{z;t) 




1 + 2^/^ 



1 -2\/< 



(5-9) 



t J \ t 

Indeed a simple Laurent expansion at cxi yields h' [z; t) — —tz^ — z+ 1 +0{z^'^), and evidently h'{^z; t) = 
^h' {z; t). Although the curve is of genus 1, the integral of h' [z) is elementary and a direct computation 
yields (recall that h{z) vanishes at one of the branch-points) 



h'{z;t)dz = h(0;t) 



yr^t 



In 



1 + 

2^/^ 



(5-10) 



We leave it to the reader to verify that dih{0; t) is continuous at i = | by using the identity ^-yi ^~^2-!/t 
1. The implicit equation of this breaking curve is then simply ^{h{0; t)) — 0. The curve is the one joining 
t = i to the t = in Figs. 4, 5, 6, 8, 9, 10. 



0' — 
-0.3 




Figure 4: All breaking curves, summarized: they are symmetric about the real i-axis. Depending on 
the case under study, some of them may not be "active" , namely they belong to different sheets of the 
phase portraits. We refer to Figs. 5, 6, 7, 8, 9, 10 for the specifics. 
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5.4 Phase portraits or distribution of genera in the complex t-plane 

There are six different situations depending on tfie values of Vj 's in the definition of the bilinear form eq. 
(1-2). Note that only their values up to common multiplication by nonzero constant is relevant, i.e., the 
orthogonal polynomials are parametrized by points [vi : vi : 1^3] e CP^. So, we have: 

1. "Generic" case: ^ 0, ^ 1^2, 1^2 7^ ^^3; 

2. "Real axis": V2 = vi, ~ 0; 

3. "Single Wedge": vi ^ Q ^ V2 ,^3 + 0. 

4. "Consecutive Wedges": 1^3 = 0, 7^ i^i, ^^2 7^ 0, v\ ^ 0; 

5. "Opposite Wedges, generic": V2 = 0, 1^3 ^ i^i, v^vj, ^ 0; 

6. "Opposite Wedges, symmetric": V2 — ^, ^ v\ ^ 0. 

We provide the results of the computer-assisted investigation for all six cases in the tables that follow 
(Figures 5, 6, 7, 8, 9, 10). The common feature is the following: as we move around the origin t = 
counterclockwise the asymptotic directions of the integration contours move clockwise by arg(i)/4. 
Therefore, a counterclockwise loop around t = yields a new configuration of contours obtained by a 
clockwise rotation of tt/2 of the initial one. 

In general, thus, we can expect that our phase portraits to have four sheets. In the "Generic" and 
"Opposite Wedges, symmetric" case, however, these four sheets are actually identical, and in the case of 
"Real Axis" two of them are equal. 

In the pictures that follow the cut (if necessary) is always along the negative i-axis and the gluing is 
the top of the negative axis of sheet j is glued to the bottom of the negative axis of sheet j + 1 (mod 4). 
We hope that the pictorial representation will serve more than many pages of verbal explanation. 

We rather explain briefly the algorithm used to investigate the phase portraits; in [2] an algorithm to 
find numerically "Boutroux curves" was explained. The algorithm produces a solution of the "modulation 
equations" (for branch-points) in high genus, but will not enforce the sign distribution (sign conditions 
for h{z)) needed to have an appropriate g-function. In terms of the Remark 4.3 it may yield a signed equi- 
librium measure. Plotting the level curves allows one to decide unambiguously whether the numerically 
produced g-function satisfies the sign distribution. 

The pictures below are produced by some code written in Matlab which is available upon request; the 
code will allow "interactive exploration" of the t-plane and to produce the pictures interactively. 

Remark 5.1 (Zeroes of the orthogonal polynomials) In all situations considered below, see Fig- 
ures 5-10, the main arcs consist of all the red arcs that are surrounded by the shaded (light blue) regions 
on both sides. These arcs, as is well known (see for example [3, 2]), also represent the limiting arcs where 
the roots of the orthogonal polynomials accumulate, and the (weak) limit of their density can be recovered 
from the jump of h'{z). 
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Figure 5: Generic. The region inside the curve joining j to is of genus 1; inside the curve that joins 
t — —-^2 ^-iid i = it is of genus (symmetric about z — Q). Everywhere else it is of genus 2, except for 
the degeneration to genus occurring on the curve that joins t = —j2 ^'^'^ ^ ~ 1' ^'^^ genus 1 on the 
ray [|;, oo). There is a Painleve I transition &i t — —j^ and a Painleve II transition a.t t — ^. 
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Figure 6: Real Axis. There are two 
sheets glued along t € K_; the level 
curves are always symmetric about 
z = 0. On the first sheet the solution 
is always of genus 0. On the second 
sheet it is of genus outside of the 
curve connecting t — —j^ and t = j. 
The region inside the curve joining | 
to is of genus 1; inside the curve 
that joins t = —j^ and i = it is of 
genus (symmetric about z = 0). In 
the remaining part it is of genus 2. 



Figure 7: Single Wedge. There 
arc four sheets glued along t G K_; 
shown here are only sheet 1 and 2, 
because the sheets 3,4 are copies of 
sheet 1, 2 where the function h{z) has 
undergone z i— > —z. Note that there 
at the critical point t = on all four 
sheets we have a transition of type 
Painleve I (and also at t — —j^)- 



Figure 8: Consecutive Wedges. 

There are four sheets glued along 
t S K_; shown here are only sheet 
1 and 2, because the remaining two 
sheets are copies of sheet 1 where 
the function h{z) has undergone z i— )■ 
—z. Note the Painleve I transition at 
both i = — and t ~ 




Figure 9: Opposite Wedges, 
generic. There are two sheets glued 
along t S K_; shown here are only 
sheet 1, because sheet 2 is a copy of 
sheet 1 where the function h{z; t) has 
undergone h(z,t) 



i-> h{z;t). Note 
the Painleve I transition a.t t — —j^ 
and Painleve II transition at t = 



1 

4- 
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Remark 5.2 Although the results, presented in Figures 5-10 are numerical, there is a straightforward 
way for their analytic justification. Consider, for example, the Generic case shown in Figures 5. According 
to Lemma 4-1, the genus zero region contains the interval {~j^,Q). Since a break can only occur at one 
of the curves defined by (5-3), the region contained inside the black curve on Figures 5 is the genus zero 
region. According to the continuation principle in the parameter space (see [3] and [19]), for any t on the 
four-sheeted Riemann surface S (with branch-points t ~ and t — oo) there exists a contour J7 = 9Jl U £ 
and hn{z), where 9Jl is the union of all the branch-cuts of hn, such that 3fJ/i„ satisfies the sign conditions 
on Q. Since there are 8 legs of zero level curves of ^th, the genus of the solution for any t cannot be 
greater than two (as there can be no bounded closed loops of^h ~ 0). Let us take a point t^, "^t^, > 0, on 
the main branch A of the breaking curve (5-3), see Subsection 5.1, that contains genus zero region inside 
(the curve from ~j2 to 0). Since is on the breaking curve, there exists a z^ £ C, such that the pairs 
{t.^,,zLz^,) satisfy (5-1) (here we use the evenness of^h{z)). Choose so that Qz* > 0. If we can show 
that 

mt{z,;t,)>0 (5-11) 

as t crosses A along 5ftt ~ 5Ri* going up, then we can prove that the genus of h is changing from zero to 
two as t crosses A. According to the Cauchy- Riemann conditions, (5-11) is equivalent to 

Qhtiz,;t,) <0, (5-12) 

where 5f = Using (4--36) and the fact that 1 -f ^ — \- t^.z'^ — (which follows from h! [z^,; t^,) — Q), 

we obtain z1 = — \ — \ , so that 




(5-13) 

To calculate the branch of the square root in (5-13) we take — > t^. As shown in subsection 5.1, in this 
limit arg(<:» — Iq) — ?• so that, using (4-16), we obtain 

111: 

arg/it(z,;i,) ^ — . (5-14) 
That proves inequality (5-12) when t^ is closed to to- Moreover, for any t^, G A we obtain 

ht{z,;t,) = ^•y/2Vl + 12i*-(l-H2t*) or /i(m) = -12f^ (5-15) 

W3tl ^ (u^ - 1)^ 

where u = \/l + 12t^ . R is easy to see that the upper halfplane part of the genus zero region (between A 
and R) is contained in the semistrip < SRu < 1, > of the u-plane. Direct calculations show that: 
'^h(u) — Q on [0, 1], and; '^h{u) < on on 1 -\- iM+ and on any segment 3m = y, € [0, 1], where 
y > is sufficiently large. Thus, using the maximum principle for Qh{u) and the fact that ^ on 
[0, 1], we conclude that 5/i(m) < inside the semistrip. So, we proved the transition from genus zero to 
genus two across A. Similar considerations will lead to rigorous proofs of transitions through other level 



33 



6 Double and multiple scaling analysis near the Painleve I gra- 
dient catastrophe points 



Having disposed of the global analysis of the problem in the complex t-plane we now focus on the so- 
called double (and multiple) scaling analysis near the two points of gradient catastrophe that are related 
to the Painleve I transcendents. These are 

io := - ^ , and — ^ ■ (6-1) 

There is another point of gradient catastrophe at t2 \ which -however- involves the Painleve II 
transcendent and should be analyzed in a separate work. 

The two points tg^ti can be analyzed much in a parallel fashion: they both necessitate of the same 
type local parametrix near one -or both- endpoints. The differences between the two cases appear by 
inspection of the figures: indeed 

• near t = to the genus zero function h{z;t) — ho{z,t) always has symmetric level-curves and hence 
the Painleve I parametrix (first introduced in [14]) is needed near both endpoints Xo — —b — — Ai 
(see Figs. 5, 6, 7, 8, 9, vignettes near to) ; 

• near t — ti the genus zero function h{z;t) = ho{z,t) does not have any special symmetry and the 
PI parametrix is needed only near one endpoint (see Figs. 5, 6, 7, 8, 9, vignettes near ti). 

6.1 Local analysis at the point of gradient catastrophe 

Near an endpoint the genus-zero /i-function has necessarily an expansion of the following form 
^hiz) = d,'\z - A,)t -f C'^\z - A,)t -f C^(z - A,)^ + • • • = 

C(^)(z-A,)t l + -i-y(z-A,) + -^(z-A,)2 + ... , J =0,1, (6-2) 

where the coefficients C^,"''' = c'^£' (t) depend on t. The gradient-catastrophe occurs when the leading 
coefficient Co\t) vanishes at one or both endpoints Ao.i of the main arc 7m, while (in general) the next 
coefficient Ci \t) does not. For our f{z,t), the gradient catastrophe point is either to or ti. Elementary 
singularity theory [1] guarantees the validity of the following definition. 

Definition 6.1 (Scaling coordinate) The scaling coordinate (^{z) — (^{z; t, N) and the exploration 
parameter t = T{t,N) are defined by 

jhiz; t) = {z\ t, N) + T(t, N)C-^ (z; i, N), (6-3) 

where C,[b] t, N) = 0, (^{z; t, N) is analytically invertible in z in a fixed small neighborhood Dj of z ^ Xj 
and T is analytic in Cq^ at Cq^ = 0, where j = 0,1. 
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Let us consider the endpoint Ai near the pomt of gradient catastrophe t^,, where = to or t^, = to. The 
expression (6-3) is the normal form of the singularity defined by h(z; x, t) (in the sense of singularity 
theory [1]). The local behaviour (we suppress the superscripts) 

2 

C = TVtQcly (^1-^^-1- 0(C2))(z-Ai)(l + O(.-Ai)), (6-4) 

r = NlCo(^^y {l + 0{Co)) (6-5) 

was calculated in [5] . The determination of the root is fixed uniquely by the requirement that the image 
of the main arc 7^, where dih = 0, be mapped to the negative real ^-axis. Following [5], we define: 

Definition 6.2 The double scaling near t — shall he defined as the appropriate dependence of t such 
that the variable 

_ 6 

v = v{t,N):=lT\t,N) = ^-mCl(^-C,^ \l + 0{Co)) (6-6) 

is kept within a disk of arbitrary but fixed (in N) radius around v — 0. The variable v shall be referred 
to as the Painleve coordinate. 



Lemma 6.1 In the double scaling near t — to for the symmetric genus zero case or near t — ti for the 
non-symmetric case, the Painleve coordinate v has the following expansion 

v{t) = ^T'{t,N)^-3hi (^t+^^NHl + 0{Vt^o)), (6-7) 

v{t)^^T^{t,N)^3i2He'^ (^t - {1 + 0{Vt^)) . (6-8) 

In either cases the function v{t) is a convergent series in y/t — tj ; if v is kept bounded as N ~> co then 
t — tj = 0{N^i). Therefore from (6-7, 6-8) it follows immediately that with accuracy 0{N^i) the map 
v(t) is linear in *-Ar . 

The proof is a direct computation with the help of Table 2 and Def. 6.1. 
6.2 Asymptotics away from the poles 

The asymptotic analysis now depends on the regions in the Painleve variable v (6-7) or (6-8) that we are 
investigating. We will split this analysis into the following two cases. 

• Away from the poles: the variable v is chosen within a fixed compact set K, that does not 
contain any pole of the relevant solutions to PI; 

• Near the poles: the variable v undergoes its own scaling limit and approaches a given pole at a 
certain rate. 
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"XT J. 1 fj- -L -L 

Near to = —j^^ ^ ^ 


M J- 1 £■-£ J- J- 

Near = ot -.^ t — ti 


a = 


a- , - 3i + Vl55t + 0((5t) 


b- -^/8 ^Vl25t + Of^t) 
Vl + Vl2(5i 2 


/)- -2i + \/l5(5t + 0((5t) 
\/l + Wlb6t 


Co ^-1(2 + ^^)725^- ^^^^^^^ 


1 

Co = — tV2ba{2a + 36) = -c^ Vl56t + 
3 3 




li(2«^ + 15a6 + 85^)^2e^^^^ 


20\/26 60v^ V / 


5 2^26 15 


N-iC{Xi) = 3^52-1^ + 0{N~i) 


N-iC{Xi) = 6"5e^ (^1 + OiN-i)^ 


CW) = + = -VsVust + 0{5t) 




3 2 
f = -- +ln2-6(5i + -(12<5i)5 +0(,5t2) 


^ ^ ^ +ln(-l) - ^(IbSt) - ^i(15<5i)' +0((5t2) 



Table 2: The explicit expressions and relevant expansions of the indicated quantities: these are the result 
of straightforward algebraic manipulations using (4-16, 4-20, 4-25, 4-28, 4-30, 4-31, 6-2, 6-3, 6-4, 6-5). 



Each of these two cases requires a slightly different analysis depending on the nature of the gradient 
catastrophe point, be it io = — or = • In. the former case the analysis was carried out in full in the 
regime "Away from the poles" by [10] and the relevant theorem is Thm. 1.1: we will not add anything 
to it. 

The only case that is not covered by the mentioned theorem in the same regime is when t undergoes 
a double scaling limit near ti and a special Painleve parametrix is needed only at one cndpoint, say, at 
Ai. Of course, one may still use the results of [10] with minor modifications to cover this new case, but 
since we will need some preparatory material, we briefly analyze this case below. We shall construct an 
approximation to the matrix T(z; t, N) appearing in (4-42) in the form 

(£(z)^'o(z) for z outside of the disks Di,Do7 
£(z)\E'o (2)7^1(2) for z inside of the disk Di, (g.g) 
£{z)-^o{z)Vo{z) for z inside of the disk Do, 
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where Dq, Di are smah fixed disks centered at Ao and Ai respectively, see Fig. 11 and as in (4-51). 
Here £{z) is the so-called error matrix that will be shown to be close to the identity matrix 1 and 'Po,i{z) 
are local parametrices at z — Ao,Ai, respectively, that will be constructed through the matrix ^'(^,t)) 
defined by (2-4). A local parametrix V{z) (we drop the indices for convenience) must have a certain 




Figure 11: The deformation of the contours and the partitioning in complementary (black) and main arcs 
(blue). Shown also are the lenses and the disks near the two branch-points Ai, Ao- The left frame refers 
to the case t ^ j^, the right frame to t ^ —j^. The thin lines in the shaded area on the left frame are 
the level-curves passing through the saddle point. For the case on the left (t = ^) we have = J^i S C, 
and i>3 can be normalized to — 1; the deformed contour 1373 consists of the complementary arc on the 
imaginary axis, the left one and the main arc (thick). The contour W2 + voi (homological sum) consists 
of the two left/right complementary arcs. For the right frame, we have t = and the weights are 

i/2 = 1 and i^i,i^3 G C: the contour vji is deformed to go through the two complementary arcs on the 
bottom and top right, and ^2 runs along the top right, top left complementary arcs and the main arc, 
while -073 runs along the two complementary arcs on the left. The weights of the various complementary 
arcs are indicated in the figure and determine the parameters of the Painleve parametrices to be used 
according to Definition 2.1. The level-curves in these pictures are numerically accurate. 



number of properties (see Theorem 6.1), one of them being the restriction 



V{z) 



l + Oe(l) 



(6-10) 
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on the boundary of the respective disk D, where o^{l) denotes some infinitesimal of e = jj, miiformly in 
z e 5D and in t e = v-'^(K). 

If the local parametrices 'Pq^i{z) satisfying (6-10) can be found then the "error matrix" £{z) is seen 
to satisfy a small-norms RHP and, thus, be uniformly close to the identity. More precisely, the matrix 
£ has jumps on: 

(a) the parts of the lenses and of the complementary arcs that lie outside of the disks Do,Di, and; 

(b) on the boundaries of the two disks Do,Di. 

The jumps in (a) are exponentially close to the identity in any norm (including L°°) while (b) on 
the boundary of the disks Dq.i we have 



= £-{l + o,{l)). (6-11) 



From the analysis in [8] it follows that, for \z\ large enough, \\£{z) — 1|| < ^^j^ (with the pointwise matrix 
norm) and that the rate of convergence is estimated as the same as the 0^(1) that appears in (6-10) as 
e ^ 0. 

In the case at hand we keep in mind that near t = ti = the endpoint Ao requires the standard Airy 
parametrix and that the corresponding error term arising on the boundary of Dg is of order 0{N^^). 

Definition 6.3 (Local parametrix away from the poles) Let C{z',£) be the local conformal coordi- 
nate near Ai introduced in Def. 6. 1 so that 

^h{z;x,t)=e{CT) = -Cl +T<:'^ . (6-12) 
2 5 

Let ^(^; v; x) denote the Psi-function of the Painleve I Problem 2.1 according to Def. 2.1. The parametrix 
V{z\ x) is defined by 



1 1 

i ~i 



r"^* (C + ^; Ir'; e-'^(f^-)'^^ , C C(^). (6-13) 
Theorem 6.1 The matrix Vi{z) := V{z]Vi) satisfies: 



1. Withinliii, the matrix Vi(z) solves the exact jump conditions on the lenses and on the complemen- 
tary arc; 

2. On the main arc (cut) 'Pi(z) satisfies 

Vi^z) = a2Vi^-{z)a2 , (6-14) 

so that within Di solves the exact jumps on all arcs contained therein (the left-multiplier in 

the jump (6- 14) cancels against the jump 0/^*0 j; 
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3. The product '^{){z)'Pi{z) (and its inverse) are -as functions of z- bounded within Di, namely the 
matrix Vi{z) cancels the growth o/^g at z = Ai; 



4-. The restriction ofVi{z) on the boundary of Di is 



= l-\Hi + —\ 



\ 0-3 1 



ley vc 2C 



H, 



16 



^2 



+ 0(r5), (6-15) 



where v — |t^, y{v) — y^^\v) and Hj = \{y')'^ + yv ~ '2y'^ ~ f y{s)ds 



The proof of Theorem 6.1 can be found in Theorem 5.1 of [5] (although it was for the tritronquee 
solution, the proof is identical for the general case).'' 

6.3 Computation of the correction near ti: proof of Theorem 1.2 
Proof of Theorem 1.2. According to (6-9), we have 

£+ = £_*o^o;i*o ^ ^£-{1 + AM{z)) on dDo^i. (6-16) 
In particular, according to (6-15), 



16 



WCp V C 2 

1- (y + ^) <J2] +0{N-i) , p:= 



z — \i 
z - Ao' 



(6-17) 



Using 1 -I- AM(z) to denote the jump-matrix of S on all the contours (see below), we can rewrite (6-16) 
as the integral equation 



£{z) = 1 



1 [ E-{s)AM{s)ds 



2i-K 



(6-18) 



where the integral is taken along all the jumps of £, that is, along the parts of the lenses and the 
complementary arcs that lie outside Dq U Di as well as along the boundaries of Dq, Di. However, the 
contribution to £ coming from the integrals along all these contours, except for 9Di, are of order not 
exceeding 0{N^^) (note that the parametrix in Dg is the standard Airy parametrix). Therefore, to obtain 
the leading order solution, we consider (6-18) with the contour dUi. This integral equation will be solved 
by iterations. The first iteration yields 

2 



£^^\z)^l + 



1 



Ai — z 



H, 



(T3 - ICTi 



1 



16/ \ 2VC'(Ai)/(Ai-Ao)/ 2C'(Ai) 



^^+i6 i-(^+D^Op 



Retaining only the terms up to order 0{N = ) in the second iteration, we obtain 



£■(2) (2) ^ S^_'\z) + res \Hi + — 



s—Xi 



.3\ ^ 



16 



(T3 - ICfi 



2VC'(Ai)/(Ai-Ao) ; (Ai - s){s - z) I y C(s) 



As 



The parametrix Vi coincides with the parametrix considered in [5] up to conjugation by cr2. 



39 



'^1 



(1) 



Ai — z 



16 



Ai — z 

(73 - iCTi 



-'^^+16 
1 



1 



1 



C'(Ai); Ai-z 

T 



1-0-2' 



2VC'(Ai)/(Ai-Ao) / 2C'(Ai) 



y 



- H 



16 



Therefore, using the fact that Xi = a + b, Xq = a — b, we have 



T(z) = 1 



Ai — z 



+ 0(iV-5) 



(Ai-Ao)a2 (Ai - Aq)^ - 4(Aj - Ag)a2 
4z 32z2 
2£'i + 6(72 (a + fe-Ei(72 6^ - 4afe(72 



2z 



2z2 



8z2 



From this we can rea<i off the relevant matrix entries: 



(Ti)22 = 



16 



V2C'(Ai)/5 
(Ti)i2=zG 



=: G, 



(Ti)2i=zG- 



2 2C'(Ai) 
2 + 2C'(Ai) 



12 



(a + 6) 



iG + 



2C'(Ai 



2/ 



ibG 



where aU the terms have accuracy 0{N 5). Direct computation using (4-48) shows 

Using Table 2, we see that 

where Oq = — 3i, 69 — 2i, and, thus 



" " 4 2C'(Ai) 



/3n — flO 



y 



C'(Ai) 



2 _3iii 

zn 5 e 10 
-1+ y{v)+0{N-^), 

5 e 10 3 
-i'i-^;^y{v) + 0{N-^). 



(72 H6-20) 



(6-21) 

(6-22) 
(6-23) 
(6-24) 
(6-25) 

(6-26) 
(6-27) 

(6-28) 
(6-29) 



To compute h„ we use (4-48). Noticing that G = 0{N 5), we can rearrange {Ti)i2 as follows: 



iba 



ly 



2C'(Ai 



iG+ —G^ + 0{N'^) 



ibo 
2 



1 



&oC'(Ai) bo 



2G 2G2 3," 



ibo 
2 



/)2 



y 

6oC'(Ai 



eT^(l -|-C'(iV-s)).(6-30) 
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Therefore 



h„ = -2z^(Ti)i2e^^ = 7rfeo(l- 



&oC'(Ai) 



e^^+T^f(l + C)(iV-t)). 



(6-31) 



Utilizing (6-22) and the values in Table 2 we find: 



2i7r 1 



2m 1 



y 

2<'(Ai) 

y 



exp 



Ne-tl + ^(C'(Ai) 



(l + 0(iV-5)) 



2iC'(Ai 



exp 



= 2i7r(-l)^^ ^1 

(l + 0(iV^5)) 



e 5 



H, 



3-i2i Ni 



exp 



V-<'(Ai) 
97V 137V 



4(15(5t)i AT _^ 



{l + 0{N-i)) 



, (15(5i)-iV^i(15(5t)5 - Jera^-^^li^- ^(155i)^ 
4 ^ ^ 3 \ ATs 3 



(6-32) 



2i7r(-l) 



N 




exp 



9iV 



195/V 



5t + 



Ni 



-H, 



(6-33) 



Q.E.D. 



7 Analysis near the poles: triple scaling limit 

The analysis in [10] was carried through under the assumption that -in the double scaling limit- the 
Painleve coordinate is chosen in an arbitrary compact set that does not contain any of the poles of the 
functions y'-^-', y^^-* (see Theorem 1.1). Our special interest now is the analysis in the vicinity of anyone 
of such poles. 

To set the stage in general terms, we shall consider the case where the Painleve variable v undergoes 
its own scaling. If Vp is the pole under scrutiny, we shall consider the following triple scaling limit, 
whereby, in addition to A'^ — >■ oo and Ni {t — t^,) being bounded, we also impose 

v-Vp^o(^N-^-py (7-1) 

where p > (depending on the situation, it may be bounded above). 

There are two distinct scenarios depending on whether the coalescence of the saddle points (zeroes of 
h'{z)) with the branch-points Aq^i occurs at both branch-points or only at one, say, at Ai. These scenarios 
corresponds to the analysis near the critical points to — ~T2^ ^^^^ *i ~ Ts respectively. We recall that 
to = is a point of gradient catastrophe in all the situations discussed in Sect. 5.4, with the exception 
of situation "Opposite Wedges, symmetric" (Fig. 10). Viceversa, the gradient catastrophe point ti = 
occurs only in "Single Wedge" (Fig. 7) and "Consecutive Wedges" (Fig. 8). 
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7.1 The asymmetric case 



Under this title we treat both the case where t is near ti (which requires a special parametrix only near 
one endpoint, say Ai, and the standard Airy parametrix near the other) and the case of t near to but 
with vi ^ I — vt, The latter case requires some special parametrix at both endpoints; but a given value of 
V, generically, can be near the pole of only one of the two special solution y^'^'>{v), of the Painleve 

I equation that enter in Theorem 1.1. Below, we assume that v is close to the pole Vp of y^^\v). The 
case when a pole Vp of y^^\v) is simultaneously a pole of t/^°'(w) even though i^i ^ 1 — t^s and, thus, 
yW ^ y^-^\ could be treated as the symmetric case (Subsection 7.2) with minor modifications (but we 
shall not consider it here for simplicity). 

We define the approximate solution to the RHP (4-2) with the jump matrix (4-43) as 



S(z)*o(^) 
E{z)^o{z)ri{z) 



for z outside of the disks Dq, 
for z inside of the disk Dq, 
for z inside of the disk Di. 



(7-2) 



where the matrix E{z), discussed below, is needed to "adjust" the situation due to the pole Vp. Here 
the parametrix Vo{z) is the Airy parametrix if we are near ti. If we are near Iq, the parametrix Vo{z) is 
given by 

Voiz) <J3Vi-z; 1 - ^.3)^x3, (7-3) 

where V{z; x) was introduced in Definition 6.3. To introduce the parametrix Vi{z), we first define ^' by 
the Masoero factorization ([17]) 



-CT3/2 



1 



*(^;i;;>f) 



(7-4) 



with 5* as in Def. 6.3 and y — y{v; x) (prime denotes derivative in v). 

Definition 7.1 (Local parametrix near the poles.) The parametrix Vi{z) is defined in . 

T 3 



Pi as 



Vi{z)=Vi{z-ui) 



1 

71 



-I 1 



(7-5) 



where C,{z\£) is the local conformal coordinate in Di, see Definition 6.3. We can then formulate the 
statement corresponding to Theorem 6.1 for the new local parametrix. 

Theorem 7.1 (Theorem 6.1 in [5]) The matrix Vi satisfies: 

1. WithinUi, the matrix Vi{z) solves the exact jump conditions on the lenses and on the complemen- 
tary arcs; 
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2. On the main arc (cut) Vi{z) satisfies 

Vi+{z) ^ a2Vi-{z)a2 , (7-6) 

so that '^qVi within Di solves the exact jumps on all arcs contained therein (the left-multiplier in 
the jump (7-6) cancels against the jump of'^^); 

3. The product "^i^){z){z — \i)^'^'^'Pi{z) (and its inverse) are -as functions of z- hounded within Di, 
namely the matrix Vi{z) cancels the growth of"^Q{z){z — Ai)""'^ at z = Xi; 



4-. The restriction ofVi{z) on the boundary o/Di is 

Vi{z) 



= (l + 0(r^)) I ^jA , (7-7) 



1 + VC/y, 



where 0(C 2 ) is uniform w.r.t. v in a small, compact neighborhood of a pole Vp that does not 
contain any zero of y{v). 

The statements in [5] were tailored to the case of the tritronquee solution and there was a slightly different 
normalization, but the proof goes through in identical fashion. Also note that the parametrix in [5] differs 
from Vi by a conjugation by 02- 

7.1.1 Triple scaling: proof of Theorem 1.3 

Before delving into the proof we make some preparatory remarks: first off, recall that we are choosing 
V so that V — Vp — 0{N~i~''), p > 0; this means that y(v) = , ^ ^2 + 0{v — v^Y ^-Iso grows at a rate 
y{v) = 0(7Vt+2p). YiecaW also that for z G 9Di we have ({z) = 0{Ni); therefore 

^^O(N-^P), zedBi,p>0. (7-8) 
y[v) 

In the case p — the disk around Ai shall be chosen sufficiently small so that \C/y\ < 1 — S for some 
(5 > 0; this means that the rightmost factor in (7-7) is a uniformly smooth and bounded matrix on 
9Di. In fact it also tends to the identity if p > 0, but in general it does so very slowly (in N) or not at 
all (if p = 0, which is the most interesting case). Therefore we can move the rightmost factor in (7-7) to 
the left at "no cost" . So, we can write 



Vi{z) 



If /9 = 0, the above mentioned factor does not tend to identity. 

We require that the approximate solution ^{z) from (7-2) satisfies 

f <I>+(z) = $_(z)(l + 0(1)) uniformly in z on aOo U 9Di 

$(7) = < (7-10) 
I $(z) is bounded for z inside the disks Dq UDi, 
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In particular, in view of point 3 in Theorem 7.1, the requirements of (7-10) will become true if the matrix 
E{z), introduced in (7-2), would satisfy the following RHP problem for E{z). 



(7-11) 



Problem 7.1 



E{z) = 0(l)(z - Ai)-'^^(cr2 + era) as z ^ Ai, 

E{z) = l + 0{1) as z oo, 

where 0(1) means an invertible matrix analytic at z — \i, hounded together with its inverse, and the 
circle dUi has positive orientation. 



Note that the second condition of (7-11) is equivalent to 



E{z)^^{z) 



-1 1 



C(z)i''" = 0(1) as^^Ai, 



(7-12) 



given that C,{z) — 0{z — Ai). Equation (7-12) together with Theorem 7.1, item 3, guarantee the bound- 
edness of 



=0{1) 



$(z) = E{z)^o{z)Vi{z) = Eiz)^o{z) 



V2 



-1 1 



=o(i) 



within the disc Di. 



Proof of solution of the Problem 7-11 Let E{z) = \{a2 + (Tj,)E{z){a2 + a^). Then 



(7-13) 



E+{z) = E^{z)M{z) on aOi, 
E{z) ^ 0{l){z - Xi)-"^ asz^Ai, 
E{z) = l + 0{-) as z cx). 



(7-14) 



where 



^(c72 + a3)^o ( *o'(^2 + a3). 

^ V 1 + V C/y / 



Using (4-51) and the fact that 



^ginSrr, = cosh(ln ^tTa) + sinh(ln fi-ffa) = ]^{S + S-^)l + ]^{S - S-^)a2, 



we calculate 



(7-15) 



(7-16) 



Miz) = [Aiz), Biz)] = , ^ ( 1 + M / — cr+ - M / — cr- I , where p - — ^ 

/i_c\ yy \i py I z- Xq 



(7-17) 
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and (7± = ^(o'l ± *cr2)- Wc make the Ansatz that E^{z) = 1 + jz^'- then the (constant m z) matrix L 
must be chosen so that E+{z) = E^{z)AI{z) satisfies 



E+{z)iz-X,r= 1 



L 



z — Ai 



[A,B](z- Ai)'"^ =0(1), zeDi. 



(7-18) 



In light of (6-4) we see that ^(Ai) = 0(1), and thus we need to consider only the second column of 
(7-18): 

L \ B{z) B{z) LB{z) 



z — Ai y z — Ai z — Ai (z — Ai)^ 

r LS(Ai) = 

1 B(Ai)+iS'(Ai) =0. 



0(1) 



Calculating 



we see that 



i?(Ai) = (0,l)^ , S'(Ai)- h. 



/ C(Ai) C'(Ai) 
2by ' 2y 



, . (Ai-Ao)y 



solves the system (7-20). Thus 



£;_(z) =l + i 



i;(Ai-Ao) 
C'(Ai) ' 



i?+(z) 



1 + i 



z — Ai 

C'(Ai) 
J2 — Ai 



A^(z) 



solves the RHP (7-14). Since |((72 -I- (T^)a^{(j2 + (73) = ^(cra — icri), we obtain that 

E-{z) = 

E+{z) = I 1+ " ^ I *n 

solves the RHP (7-11). Q.E.D. 



2(z-Ai) 
2(z-Ai) 



(7-19) 
(7-20) 

(7-21) 
(7-22) 



(7-23) 



(7-24) 



Error analysis. The error matrix 8{z) — T(z)$~^(z) has jumps on the lenses and on the complemen- 
tary arcs outside the disks Do,Di, as well as on the boundary of these disks. The jump matrices on the 
lenses and on the complementary arcs approach 1 exponentially fast in N^^ and uniformly in z. It is 
also clear that 5 — > 1 as z ^ 00 since both T and $ = E'^E'o do so. So, it remains only to prove the 
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uniform convergence to 1 of the jump matrix on 9Di (convergence on 5Do was established in Subsection 
6.2). Indeed, using (7-2), (7-11), (7-7), (7-9) and (7-24), we have 

-0-3 



,1 + VC/y, 

f_£;_«'o(l + 0{C^''^))-^^^EZ^ = £-E_{l + 0{C,-^''^))EZ^. 



(7-25) 



On the boundary z G 9Di we have C, = OiNs) and in our triple scaling y — 0{N5+'^p) with p > 0. Then 
L is of the order 0{N-t,^) = 0{NP). Thus, 

z-Ai [z-Xiy 



So, it is the last term that contributes the slowest decay. Therefore, we obtain 
The latter estimate shows that we can control the error provided 



y = = 0{Nt+P), or, equivalently, v - Vp ^ 0{N~'s-p), 



(7-27) 



(7-28) 



where < p < i . 



Computation of the recurrence coefficients: We need to use (3-9) and (7-27). Using (7-2), (7-24) 
and the expansion of '^q (4-52) we obtain 



Hz)^E^iz)Mz)={l+l'-^^^^^Ul 



2 z-Xi 



b abu2 ^^ 

^-2 + ^l-^+0(.-3) 



z z^ 



1 (To 

2z ^ 



fe^ aba-;, o 



= 1 



k := 



1 



- (cTg - ICTl) - -(72 



1 



2z2 

b^l ab k f 1 « , 

^"y"^ + 2 1"+2^'^"3-^"^) 



/(Ai-Ao)y _ 26y 



C'(Ai) 



C'(Ai) 



We introduce 



s 



(7-29) 
(7-30) 

(7-31) 
(7-32) 



where the latter expression follows from (1-14). Here and henceforth ao,&o denote the values of a, 6 
calculated exactly at one of the critical points <o or ti. Assuming in p = in (7-28), we obtain s = 0(1) 
as TV —> 00. On the other hand, if p € (0, |), then, consequently, s scales as 0{N~p). Thus, we have the 
triple scaling limit 



t-ti = 



2bo 
C'(Ai) 



(7-33) 
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where k is the constant in front oi St = t — tj appearing in formulae (6-7, 6-8). Exphcitly, using Table 2, 
we obtain 



t 



1 Vp s 

12 ""sMivt "ayfiv ' 

Now, according to (4-48), (4-30), (4-31) and (7-27), we obtain: 



1 VpQ s . 2s 

15 " 3t2s57Vt ~*15]V' 



1 

4^ 



2 +0(iV-5y,7V-5), 



1 



-2i7r(ri)i2e 



Nl 



2s(l - hos) + 0{N-iy) 
( 5o - ^ ) +0{N~'iy,N-i 



+ 0{N--2y,N-t), 



(7-34) 

(7-35) 
(7-36) 
(7-37) 



Here ©(iV-s) error term comes from replacing a,b with their respective values ao,6o considered at the 
critical point to or ti. Note, however, that in the regime (7-1), the 0{N^i) term is of a smaller order 
than the 0{N~iy) term. Therefore, in all these expressions, in the regime (7-1), the error is at best 
0{N^^) (recall that y — 0{Ni^'^P) p G [0, ^)). Thus in the exponent e^^ we can use the expansion in 
Table 2 up to order St included. So, 

N 



exp 



-iV- 



cNSt 



bo 



1 



+ OiN-^y) , 



(7-38) 



where c — —6 for the case t to and c 
expressions in (7-34). So, in the leading order. 



1315 



for the case i ~ ti. One has then to replace 5t by the 



h„ = TT ( V8 - ^ ) 2^ cxp 



TT 2i - - (-1)'^ exp 



3A^ N^Vp 



\/2i 



t 



12' 



97V 13 7V5Upe- 



3525 



13 



t 



1 

15" 



(7-39) 
(7-40) 



It is remarkable to note that the genus zero leading order asymptotics a„(i) ~ ^ and /3„(i) a are 
valid as long y = o{Ni) with the accuracy 0(-^). However, when y — 0{N^), both terms in (7-35), 
(7-36), contribute to the leading order, whereas, when y = 0{Ns+'^p) with p e [0, ^), the asymptotics are 
determined by the latter terms of (7-35), (7-36). In this case, both a„ and /3„ are unbounded as — >■ oo. 

So, the proof of Theorem 1.3 is completed. Q.E.D. 



7.2 The symmetric case: proof of Theorem 1.4 



12 ' 



We are now in the symmetric situation and hence the critical point to consider can only be to 
where Ai = 6, Aq = —b and a = 0. This case is significantly different from the previous inasmuch as the 
two Painleve parametrices in Dq^i are identical: in particular y*^^^ = y^'^-' = y. Thus, if the double scaling 
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is such that we are close to a pole Vp of y{v), this will simultaneously affect the both parametrices and, 
as we shall see, will have a significant effect on the asymptotics of «„. On the other hand, due to the 
exact symmetry of the bilinear pairing, the orthogonal polynomials have the same parity of their degree 
and thus automatically /?„ = 0. 

It will be advantageous for us to use a different solution to the model problem (4-50), which has a 
different growth rate near the branch-points: such modification (see [11]) is called a discrete Schlesinger 
transformation. In terms of the RHP (4-50), this amounts to replacing the solution (4-51) with 

^i(.):=2(^3 + ..)(^) (-3 + ..)^(^) . (7-41) 

This matrix satisfies all the conditions of the RHP (4-50) except the last one, as it clearly has a different 
growth behaviour near the endpoints ±6. We then shall construct an approximate solution 

E{z)'^i{z) for z outside of the disks Do,Di 
<i>(z) = I E{z)^i{z)'Pq{z) for z inside of the disk Do, (7-42) 
E(z)^i{z)Vi{z) for z inside of the disk Di, 

where Vi{z) is defined by (7-5) and 

Vo{z)^(T^Vi{-z)a:, . (7-43) 

Due to the fact that we are using 'I'l instead of 'So, the boundedness of the product at Ao,Ai 

follows immediately (see also Theorem 7.1, item 3). Hence, the requirements on the left multiplier E{z) 
are now different compare with the asymmetric case studied above (we reuse the same symbol E with a 
new meaning relative to the previous section). 

Problem 7.2 Find the matrix E{z) is analytic (together with its inverse) on C\(9DoU9Di) and satisfies 

E{z)^l + 0{1) as z^ 00, 

where the contours 9Do,i have positive orientation. 



E+{z)=E^{z)^,{z)[^^-^] ^Sj^\z) ondDo.i, 



Proof of solution of Problem 7.2. Note that any solution to this RHP has unit determinant and 
hence its inverse is also analytic and bounded. As before, we find it more convenient to solve the RHP 
for E{z) = FE{z)F instead of the RHP (7-44). Here F = The jump matrix for the new RHP is 

M{z) = F-Hi{z)Q-\z)-^^\z)F, (7-45) 



where 
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y = y{v), V was defined by (6-7) and the local scaling coordinate C — ({z, N) near z = b was introduced 
in (6-4). 

Direct calculations yield 



M{z) = 



1 



Similarly, near z = —b we obtain 



i-az)/y 



1 

C(^)(^-fc)^ 



1 

C(z){z-b)3 
y{z+b)3 



C{z){z+b)3 
y(z-bY 

1 



y(z-br 

1 



(7-47) 



where C = C(^) ~ C(^^) i e 



(7-48) 



(7-49) 



Note that the orthogonal polynomials in this case are even/odd and the symmetry of the RHP implies 
(which can be verified directly from the above formulae and also as a consequence of (7-45)) 



^-1(2) ^ CT3^'l(-2)cr3 

Using (7-47), (7-48), (6-4), (7-49), we obtain 



M{z) = a2M{-z)a2 



M{z) 



W+ 
z-b 



+ Oi{z), M{z)-1 



z + b 



+ Oo{z), 



where 



77 = iN ^ 



'{2bfC 



(7-50) 



(7-51) 



(7-52) 



Here C = N^^C'ib) and C = 7V"s^'(_6) = -C* by the symmetry of the problem. Note that the Oo,i 
terms in (7-51) are analytic at z = Ao,i and when evaluated at z = Ao,i = ±6 are proportional to a± 
(respectively). The matrix E{z) satisfies 



m = 

We pose the Ansatz 



E^{s){M{s) - 1) ds 



\s — b\—r 



2i7r 



S_(s)(M(s) - 1) ds 



E_iz)^l 



\s+b\=r 

A A 



2iTT 



z — b z + b^ 
and obtain 

A A _ Ar]a+ T]a+ Aria+ Arja- AOi{b) AOoi-b) -qu- A-qa. 



z — b z + b 



{z-bY z-b 2b{z-b) (z + 6)2 



-b 



z + b z + b 2b{z + b) 



(7-53) 



(7-54) 



(7-55) 



That leads to the following system for the unknown A^ A (recall that Oi{b) cx a+, Oo{~b) cx cr_): 

A(J+ = 0, Acr_ = 0, 
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A 



^ A 



A 



2b 



This system has the solution 
A ' 



1 



(2b)2 



-T] 





A 



1 



(26)2 



Aa^ = — ?7(T_ 



^ 
rj 



(7-56) 



(7-57) 



So, we found £^(2) and, thus, E{z). Note that the function E{z) in the region outside of the disks is 
a rational function with poles at ±6, while, inside the disks, it is analytic and given by formula (7-53). 
Q.E.D. 



Error analysis. The error matrix £{z) — T(z)$^^(z) has jumps on the lenses and on the complemen- 
tary arcs outside the disks Do,Di, as well as on the boundary of these disks. The jump matrices on the 
lenses and on the complementary arcs approach 1 exponentially fast in N and uniformly in z. It is also 
clear that £ — > 1 as z — > cx3. So, it remains only to prove the uniform convergence to 1 of the jump 
matrix on SDq.i: the computations are absolutely parallel and we report only the one for dDi. Using 
(7-2), the solution to Problem 7.2 and eq. (7-7), we have 



,1 + VC/y, 



-^i^EZ^ 



(7-58) 



On the boundary z e dDi we have C = 0{Nt) and in our double scaling y — 0{N^'^'^p), where p e [0, |). 
Moreover, E^ = 1 + ^jz^ + ^j^, where A, A are of the same order. That creates the situation that is 
drastically different from the previous: for example, the matrices A, A (7-57) remain bounded no matter 
how fast y grows (and hence — > (7-52)). The only unboundedness occurs when the denominators in 
(7-57) vanish, which means that rj has a finite value 77^ — —46^ or, equivalently. 



iV-5y(„) = 2Cb. 



(7-59) 



Condition (7-59) identifies two points near the pole v = Vp at a, distance of order 0{N b). Thus, in 
(7-58) we have 



E_0{C^'^)EZ^ = 0{C 



-1/2N 



462 



462 



(7-60) 



The very last contribution to the error term comes from the denominators of the matrices A^ A (7-57) 
and prevents us from getting close "too fast" to the points where they vanish. 
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Computation of the recurrence coefficients: Following [5], we find the expansion of the matrix 

^(z) = E{z)^'i{z) at z = cxi: 



E^i = FE{z) 

Using (7-57), we obtain 

A + A = - 



F ^ = 1 + ^ ? + ^-^ V 

z 



(73 + 770-1 



A-A = - 



1 - ir]a2 



+ 0{z-'). 

(7-61) 

(7-62) 
(7-63) 



so that 



F{A + A)F-^ = 



1 _i_ "1 

^ (26)2 



2b 



(72 + Wl 



F{A-A)F-'^ = 



.(7-64) 



It follows from (7-61) and (7-64) that the residue of $ at infinity, which we denote by <i>i, is 



$1 = F{A + A)F-' + h<j2 = b^^ 

^ + 462 



4fe2 



1 

1+^ 

^ + 462 



-02 



(7-65) 



We note in passing that $1 is ofF-diagonal and $2 is diagonal (which implies /3„ = 0, which -of course- is 
identity and not just an approximation due to the special symmetry of this case). Then 



($1)12 = 



($1)2 



3i + 


462 + ^ 6 _ 


63 + 

—i 


ir\ 
26 


1 - 


1 

^ (26)2 


21- 


ir\ 
26 


bii 


iri^ 19'; 
462 -1- ^ 


53- 
= i 


ir\ 
26 


2 


^ ^ (26)2 


21 + 


H) 

26 



(7-66) 
(7-67) 



Using (7-60), we can now calculate the (leading order) final expressions 

«n = (ri)i2(ri)21 = 



7,2 q 1 iL 

^ + 462 



(7-68) 



where it is understood that both expressions (also in the denominators) are affected by an error of the 
order indicated in (7-60). Introducing 



11] 
'2b 



y 



(7-69) 
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we note that 0{N-r,{l + 7j^/{4P))-^) = 0{N-i{s^ - and we find finally (using Table 2 for the 

symmetric case)^ 



4 1 - s2 + 0{N 
h„ = 2^7r\/8exp 



(3n = 0, 



3^ ^,.1 St 
— -6iVs 

2 N-t 



Using (6-7) and Table 2 to relate s and t, we can write (7-70) as 

'3 - s 



h„ = 7r%/8 2^exp 



3N Nr, Vp s 



1 + s 



O 



1-S2 



1-S2 



(7-70) 



(7-71) 



Q.E.D. 
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